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Abstract 

We investigate dynamics near Turing patterns in reaction-diffusion systems posed on the reai 
line. Linear analysis predicts diffusive decay of small perturbations. We construct a "normal form" 
coordinate system near such Turing patterns which exhibits an approximate discrete conservation 
law. The key ingredients to the normal form is a conjugation of the reaction-diffusion system on the 
real line to a lattice dynamical system. At each lattice site, we decompose perturbations into neutral 
phase shifts and normal decaying components. As an application of our normal form construction, 
we prove nonlinear stability of Turing patterns with respect to perturbations that are small in 
L 1 n L°°, with sharp rates, recovering and slightly improving on results in [21, 10]. 
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1 Introduction 



Turing predicted that the simple interplay of reaction and diffusion can lead to stable, spatially periodic 
patterns [22]. His ideas proved quite influential in the general area of pattern formation, where one 
seeks to understand the formation and dynamics of self-organized spatio-temporal structures. One 
can easily envision simple reaction-diffusion systems with two species that exhibit diffusion-driven 
instabilities of spatially homogeneous equilibria. Typical examples are activator-inhibitor systems such 
as the Gray-Scott or the Gierer-Meinhard equation; see for instance [15, 16]. Perturbations of the 
homogeneous unstable equilibrium grow exponentially at an initial stage, with fastest growth for distinct 
spatial wavenumbers. This wavenumber is roughly independent of boundary conditions in large enough 
domains. As a final result, one often finds a spatially periodic pattern, up to narrow, exponentially 
localized boundary layers. In order to understand such nonlinear spatially periodic patterns and the 
process of wavenumber selection, one is therefore naturally led to considering reaction-diffusion systems 
on idealized unbouned domains. 

To fix ideas, consider 

u< = Mu + f(u), 

for u(t,x) £ M n , with x £ WL N , with smooth reaction-kinetics f and positive diagonal diffusion matrix 
D = diag (dj) > 0. Here, and in the following, the term "smooth" refers to functions with sufficiently 
many derivatives. In many circumstances, one can show that there exist families of spatially periodic 
striped solutions, 

u(t, x) = u*(foci; k), u*(£; k) = u*(£ + 2tt; k), 

parameterized by the spatial wavenumber k > 0. In fact, such families occur for an open class of 
reaction-diffusion systems, including but not limited to systems of activator-inhibitor type mentioned 
above. 

As a first predictor on the stability of such solutions with respect to perturbations, one analyzes the 
linearization, 

v t =DAv + f'(u i (b;i))v. (1.1) 

It turns out that, again for open classes of reaction-diffusion systems including the above examples, 
solutions to this linear equation are bounded for bounded initial data, for an open subset of patterns 
u*(-;fc) in the family. We refer to such patterns as linearly stable Turing patterns. We will discuss 
detailed assumptions that guarantee such linear stability later in this section. 

The presence of a family of patterns, parameterized by the wavenumber, and, even more obviously, by 
translations of the pattern in x, implies that solutions to (1.1) with general initial conditions will not 
decay. More explicitly, v(t,x) = d x u±(kx;k) and v(t, x) = ^rU*(fcx; &;) are constant in time and solve 
(1.1). 

In fact, one can show that under typical assumptions, initial conditions v(t = 0, x) £ L 1 (R Ar ,M n ) 
will give rise to diffusive decay, sup^ |v(i,a;)| ^ Ct~ N / 2 . Such algebraic decay is in general not strong 
enough to ensure nonlinear decay in dimensions iV ^ 3. The simplest example is the nonlinear heat 
equation 

ut = Aw + u 2 , 

which exhibits blowup of arbitrarily small, smooth, positive initial data at finite time in dimensions N ^ 
3 [6, 3]. In the seminal paper [21], Schneider recognized that diffusive decay near Turing patterns is not 
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altered by the presence of nonlinear terms due to cancellations in a Bloch-wave expansion. He studied 
the most difficult case, N = 1, where diffusion is weak and nonlinearity potentially most dangerous, in 
the specific example of the Swift-Hohenberg equation. His proof has later been generalized, simplified, 
and adapted; see [23, 7, 8, 11, 9, 10, 4, 19]. Our focus here is, again, on the one-dimensional case, in 
a general reaction-diffusion setting. Our goal is to find coordinates that show explicitly why nonlinear 
terms do not alter linear decay near Turing patterns. Going back to the scalar heat equation, the 
interaction of nonlinear terms with diffusion can be categorized as relevant, critical, or irrelevant; 
[1, 2]. Explicitly, in the heat equation ut = u xx + f(u,u x ,u xx ), 

(i) Nonlinear terms such as f(u,u x ,u xx ) = uu xx ,u 2 ,u p , where p > 3 are irrelevant; 

(ii) Nonlinear terms such as f(u, Ux j Uxx ) — UUx j U^ 3X6 critical; 
(iii) Nonlinear terms such as f(u) = u 2 are relevant. 

Without pretending to fully explain this phenomenon, notice that, for LMnitial data, assuming Gaus- 
sian decay, we find u xx ~ i _3//2 in L°°. Irrelevant nonlinear terms decay with rate t~ a , a > 3/2, critical 
terms have a = 3/2, and relevant terms have a < 3/2. 

Perturbations v of Turing patterns solve a system 

vt = d xx v + f'(u*0r))v + g(x, v), 

where g(x,v) = 0(|v| 2 ). Note that from here on, we fix the wavenumber k = 1, without loss of 
generality, and write u*(x) := u*(x;l). In particular, the nonlinearity g has potentially dangerous 
quadratic terms. Roughly speaking, our goal is to find coordinates in which the nonlinearity involves 
at least two "derivatives", which according to the numerology for the scalar heat equation would 
be sufficient to guarantee nonlinear decay. The reason to hope for derivatives is the presence of a 
conservation law associated with the translation symmetry, which in turn generates the neutral decay 
in the linearization. 

To be precise, we now consider reaction diffusion systems 



where u, f £ M n ,a; 6 M,t 6 (0,+oo),D G ]R rixn is a diagonal matrix with strictly positive diagonal 
entries and f is smooth. Firstly, we assume the existence of a Turing pattern of the system. 

Hypothesis 1.1 (existence) The system of ordinary differential equations Dd xx u+i (u) = possesses 
a smooth periodic even solution u*. 

Without loss of generality, we assume that the period is 2tt. Our aim is to study nonlinear stability of 
this temporal equilibrium under general small non-periodic perturbations. To this end, we introduce 
an initial condition 



Then assuming that u(t, x) = u*(x) + v(t, x) is a solution to (1.2) with the given initial condition (1.3), 
we have 




(1.2) 




(1.3) 




(1.4) 
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(1.5) 



where 

A: X 1 — ► X 

v .— ► Da«v + f'K)v. 

Here we define 

X = (L 1 (R)) n n(L°°(R)) n , X 1 = (W r2 ' 1 (R)) n n (VF 2 '°°(IR)) n , (1.6) 

with norms 

II ' ll-X" = II ' Hi 1 + II ' II ' Wx 1 = II ' llw 2 . 1 + II ' ||vK 2 '°°- 

Note that from now on, we suppress n and M. if there is no ambiguity. Moreover, g : x W l —> W 1 is 
smooth, g(x, v) = f (u* + v(x)) — f (u*) — f'(u*)v, so that g(x, 0) = and d v g(x, 0) = 0. 

According to Bloch wave decomposition, let us introduce the family of Bloch operators, for a G [— i, g], 

B{a): (tf 2 (T 27r ))« — > (L 2 ^))" 

L>(^+ia) 2 v + f'K)v. 1 ' 1 



For further reading on Bloch wave decomposition and Bloch operators, we refer to Section 6.2 and [17]. 
Note that one obtains B(o~) formally by applying A to functions of the form u = e iax v. 

Hypothesis 1.2 (spectral stability) The family of Bloch wave operators B{a) has the following 
properties. 

(i) spec(£(o-)) f|{#eA ^ 0} = 0, for a / 0; 

(ii) spec(i?(0)) P|{i?eA ^ 0} = {0} and is simple with spanju^} as its eigenspace; 

(Hi) Near a = 0, the only eigenvalue X is a smooth function of a and the expression of \(o~) reads: 
A(cr) = —da 2 + 0(|cr| 3 ), where d > is a constant. 

Remark 1.3 The expansion in (Hi) is a consequence of the simplicity of A = at a = and the 
evenness of u+. In fact, we have an "explicit" expression for d; see Section 6.4. 

Given the above hypotheses, we can state our main result. 

Theorem 1 (nonlinear stability) Assume Hpotheses 1.1 and 1.2 hold. There are C,a > so that, 
for any ||v°||x < <? , where X = (L 1 (R)) n n (L°°(R)) n , the solution v(t) to the system (1.4) exists for 
time t £ [0, 00) and satisfies the estimate 

l|v(t)ll(L~(R))" ^C-l^fr. (1.8) 

(l+t)3 

The rest of the paper contains three main contributions. First, we construct normal form coordinates, 
where the neutral mode is represented by a discrete phase 6j , which decays according to a linear discrete 
diffusion equation 6j = d{6j+\ — 26j + 0j-i)- The idea is to capture the leading order dynamics of 
perturbations using an ansatz of the type u(i, x) = u+(x — 6j) + Wj(t,x) on intervals x £ [27r(j — 
1/2), 27r(j + 1/2)], where Wj(i, x) lies in a linear strong stable fiber. The coordinate change mimics the 
much simpler coordinate change in [4], where strong stable fibers of a temporally periodic, but spatially 
homogeneous solution were straightened out. 
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Our second main contribution are decay estimates for the linearization in these coordinates. In partic- 
ular, we show that the Wj indeed decay with higher algebraic rate than the 6j. 

Our third main contribution is the computation of nonlinear terms in the new coordinate systems. 
Leading nonlinear terms turn out to involve discrete derivatives, associated with the discrete transla- 
tional symmetry near the periodic pattern. Similarly to the scalar case, these discrete derivatives render 
the nonlinearity irrelevant. Prom a different view point, dependence on derivatives, only, indicates the 
presence of a conservation law: An equation ut = u xx + f(u x ) can be rewritten as vt = v xx + (f(v)) x , 
for v = u x , and the gain in decay is now clear from an integration by parts in the variation of constant 
formula. An analogous observation applies to the 9 — W system, where discrete derivatives in the 
nonlinearity reflect a discrete conservation law. 

Together, these observations quite readily imply a nonlinear stability result-Theorem 1 as shown above. 

The remainder of this paper is organized as follows. In Section 2, we construct the normal form. 
Section 3 contains linear estimates in Fourier-Bloch space. Section 4 converts those decay estimates 
into LP — L q decay estimates in physical space. Section 5 contains the proof of the nonlinear stability 
result. We relegate a detailed description of the nonlinearity, and the spectral properties and the 
analytic semigroup results of the linear operator to the appendix. 

Notation Throughout we will use the following notation. 

• (•, •) is the standard inner product on W 1 given by 

n 

(u, v) = ^2u jVj , for any u = { Uj }] =1 ,v = {vj}] =1 £ R n . 

3=1 

• (■, •) is the standard inner product on the Hilbert space (L 2 (— n, vr)) n given by 

(vl(x), v(x))dx, for any u, v G (L 2 (— ir, vr)) n . 

• ((•, •)) is the standard inner products on (£ 2 ) n , or the (£ p ) n -(£ q ) n pairing, given by 

fav» = ^(u i; Vj), for any u = {uj} jeZ , v = {fjjjez- 

We denote the Euclidean norm in Euclidean spaces as | • |, the norm in a general Banach space X as 
|| • \xi an d the norm of a linear operator from a Banach space 3C to as HHH^-^^- For the case 
W = JT, the last norm notation simply becomes ||| - |||^r. 

2 Normal form 

As we pointed out in the introduction, the linear system for the perturbation 

v t = d xx v + f'(u*(x))v, (2.1) 
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is expected to exhibit diffusive decay for the linear part. This weak decay is not obviously strong enough 
to conclude nonlinear decay because of quadratic and cubic terms in the nonlinearity. Our approach 
here converts the system (2.1) into an infinite-dimensional lattice dynamical system for V = {Vj}j & %, 
where V,- = (Oj, W,) Six L p (—tt, tt) for all j G 7L. Here, the scalar component Oj of Vj measures local 
shifts of the primary periodic pattern, and the infinite-dimensional component, W,-, represents local 
complements. In such a representation, one expects diffusive decay of Oj and faster decay of Wj. We 
will make this precise in Section 4. In fact, the linear asymptotics of Oj are equivalent to the discrete 
diffusion 

9j = d(9j+i — 26 j + Oj—i). 

The key idea is that in this lattice system, nonlinear terms in the ^-equations involve discrete derivatives, 
8j + i — 8j, rather than Oj alone. Roughly speaking, we expect ^-dependence to disappear when Oj = Qj+i 
for all j E Z due to shift invariance of the original system. Just like in the continuous scalar heat 
equation, these derivatives decay faster, so that terms like {Oj+i — 9j) 2 are now irrelevant, that is, they 
do not alter linear diffusive decay. 

In summary, we will find a system, where the linear part exhibits diffusive decay, and where nonlin- 
earities are explicitly irrelevant. In this sense, our transformation has eliminated lower-order terms in 
the system, that turn out not to contribute to leading order dynamics. The term normal form alludes 
to this elimination of lower-order terms by comparing with normal form theory in ODEs, where coor- 
dinate changes are used to simplify equations and systems at least locally, mostly through removing 
lower-order terms in the Taylor jet of the equation or system. 

The remainder of this section is organized as follows. We discuss local well-posedness and "chopping- 
up", the first key step in the transformation to a lattice system in Section 2.1. The ultimate transfor- 
mation towards a quasilinear lattice dynamical system is constructed in Section 2.2. Key steps involve 
separation of the neutral phase Oj and a smoothing procedure at the chopping boundaries. 

2.1 Well-posedness: spatially extended system and lattice system 

We first show local-in-time well-posedness of the system (1.4) on the space X = L 1 n L°°. 

Lemma 2.1 The initial value problem of the semi-linear parabolic system (1.4) is locally well-posed in 
X . To be precise, the following assertions hold: 

(i) (existence and uniqueness) For any given v° E X, there exists some T > 0, depending only 
on 1 1 v° 1 1 jsc, such that the system (1.4) admits a unique mild solution 

v g c°([o, t], (L\R)) n ) n c°((o, r], (L°°(M)) n ). 

Here a mild solution solves the integral- equation variant of (1.4). 

(ii) (regularity) The solution v(t,x) to (1.4) is smooth for t G (0, T\. Moreover, there exists C > 
such that, for all t G (0, T], 

\\v(t)\\ H 2 ^C^Hvlx. 

Proof. The existence and uniqueness follow directly from [5] and [12]. To show that ||v(£)||#2 ^ 
Ct" 1 1| v° || j\: , we first note that for any To G (0,T), there exists C(Tq) > such that 

||v(t)|| H2 ^ C(T )||v°|| L2 , for all t G (T ,T). (2.2) 
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Moreover, by [12, Thm. 7.1.5], for every v° G L 2 , there are Ti > and C(Ti) > such that 

\\v(t/2)\\ m ^CiT^t/ir^W^Wv, \\v(t)\\ H2 ^CiTiXtffl-WWvQMWHi, faraUte(0,Ti), 
which implies that 

l|v(i)ll^ < ^f^ _1 ||v || L2 , for all t G (0,7i). (2.3) 
Combining (2.2) and (2.3), we conclude our proof. ■ 

Remark 2.2 5y Lemma 2.1, we can assume without loss of generality that, in the proof of Theorem 
1, the initial perturbation is small in X n H 2 . 

Now suppose that v(i, x) is a solution to (1.4), close to 0. In particular, v(t, x) is close to on all 
intervals [2ir(j — 1/2), 2ir(j + 1/2)], j G Z. Then instead of solving (1.4), we claim that it is equivalent 
to solve the infinite-dimensional system, for all j G Z, 

d t \j = Dd xx Vj + f'(u*)vj + g(x, Vj) 

v i (t,7r)=v i+ i(t,-7r) (2.4) 
d x Vj(t,n) = d x v j+ i(t,-x). 

In order to justify the well-posedness of (2.4), we first introduce the chopped space 

X ch = e\Z, (L^-vr, 7r)) n ) n r°(Z, (L°°(-7r, tt)) b ), (2.5) 

with the norm defined as 

l|w||x ch = ^ II w j Hz, 1 + SU P ll w ilk°°) for an y w = {wj}jez g x ch . 

We then consider the chopping map 



-^ch '■ X c h > X 

X 1— > ^ch(v), 



A ch * A ch, 



where := fT^X 1 ) and 



(2.6) 



where X is defined in (1.6) and =^h(v)(27rj + x) = Vj(x),for all x G [— vr,7r] and j G Z. It is not hard 
to see that is an isomorphism and thus we have the diagram 

x 1 A x 
<%h t t 



^ch : > -^ch 

V .— > ^^^chX. 

More specifically, (A;hX).j = Dd xx Vj + f (u*)vj. To describe X^, we define 

®(A ch , X ch ) := e^W^i-Tr, vr)) n £°°(W 2 >°°(-tv, vr)), 
^(^l ch ,X ch ) := {v G ®(A ch ,X ch ) | vf \t,7r) = vf^(t, -n),t > 0, j G = 0, 1}. 



(2.7) 
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Lemma 2.3 We have X\ h = @(A ch ,X ch ). 

Proof. From the definition, we find X^ h C @(A c h, X c h). We only need to show that for any given 
v G St(A c h, Xck), we have v = <%h{y.) G X 1 . In fact, for arbitrary w G C£°, we obtain 

(v,w') L 2 (R) = ^(vj(x),w'(27rj + x)) = - ^(Vj(x), w(2vrj + x)) = -(^ c h({vj}j G z), w) L2(R) , 
iez jez 

which shows that v' = =%h({ v j}jez) G X. Similarly, we have v" a = =%h({ v j Ijez) G X. ■ 

In all, we conclude that our initial value problem for a spatially extended system (1.4) is equivalent to 
an initial value problem for a lattice system as follows. 

dtY. = -4chV + G(v), x G (— 7r,7r),t > 0, 
< vf\t,ir) = vf^t, -vr), k = 0, l,j G Z,t > 0, (2.8) 
>Vj (0,x) = v°(27rj + x), x G [-7r,vr],i G Z, 

where G(v) = {g(x, Vj)}j £Z . 

Remark 2.4 For any solution v to (2.8), we have all higher matching boundary conditions, that is, 
d™Vj(t, 7r) = 8J*v i+ i(t, -vr), /or a// 1 > 0, j G Z and m G Z + . 



2.2 Lattice system: phase decomposition and boundary-condition matching 

We start with sketching the construction of the normal form step by step without rigorous justification. 
We first decompose each 2-7r-long piece Vj(x) = v(27rj + x) into a linearly neutral phase and a stable 
phase and then match the boundary conditions for the stable phase. This two-step smooth phase 
decomposition procedure will be summarized and justified rigorously in a lemma at the end of this 
section. 

We now decompose each v,- according to 

jvj(x) = Wj(x) + u*(x - Oj) - u*(x) 
[(wj(x),u ad (x - Oj)) = 0, 

where u ad is an element in the kernel of the adjoint operator of -B(O) with (u^,u ad ) = 1. Substituting 
this expression into (1.4), we can therefore formally derive a system for 0j and Wj, which takes the 
explicit form 

• 1 



-1 + (W,(x),u' ad (^)) [ - (W ^ (7r) - W i(--)^ U 'a> " m 

+ (d x Wj(ir) - d x Wj(-ir),Du ad (Tr - 0j)) + (g(0j, Wj), u ad (x - 0j))\ 
k Wj =Dd xx Wj + u* >e (x - 0j)0j + f(wj + u*(x - 0j)) - f(u*(x - 0j)), 
with boundary conditions 

f ^Wj(Tr) - a> i+1 (--7r) = ui m) (7r - j+1 ) - ui m) (tt - 0j), for m = 0, 1 
1 (Wj(x),u ad (a; - 0j)) = 0, 

where 

g(0j, Wj) = f(wj + u*(x - 0j)) - f(u*(x - 0j)) - f (u*(x - 0j))wj. 



(2.9) 
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Remark 2.5 In the second equation of (2.9), 9j represents the right hand side of the first equation. 



Note that Wj is in a codimension-one subspace depending on 9j. More formally, we mapped every Vj 
into a vector bundle. Also, the boundary conditions are now nonlinear. These facts generate technical 
difficulties so that we find it easier to work with a further modified system, where, for all j € Z, we 
substitute 

Wj (x) =W,(x)+H(x,9 j - 1 ,e j> 0j +1 ,W j ). (2.10) 

For simplicity, we denote Hj(x) = H(x, Oj—i, 9j, 9j+i, Wj). In the new coordinates V = (9, W), where 
9 = {9j}j£z and W = {Wj}j g z, we will have again "homogeneous matching boundary conditions" 
and all W,-'s are in a fixed codimension-1 subspace, that is, for all j G Z, 



dgWjfr) - d™W j+1 (-n) = 0, for m = 0, 1, 
(W j (x),u ad (x)) = 0. 



(2.11) 
(2.12) 



We now construct H = {Hj(x)}j 6 ^ explicitly in the form 



H 



(2.13) 



where Hj accomplishes "homogeneous matching boundary conditions" (2.11) and H| corrects so that 
every Wj is perpendicular to u a( j (2.12). First, we construct Hj. In order to accomplish (2.11), one 
readily verifies that we need 

<5™H)(tt) - d™H} +1 (-vr) = ui m) (^ - 9 j+1 ) - u£ m) (vr - for m = 0, 1, 

which can be achieved by choosing 

H](x) = ^(u*(x - 9j+i) - u*(x - 9j)), for x ~ it, 
H}(x) = ^(u*(x - 6»j„i) - u*(x - 9j)), for x ~ -it. 



In light of this observation, we let 
Hj(x) = -0(x)(u*(a; - fy+i) ~ u *( x ~ fy-O) + t( u *( 3 



2u*(x-%)), (2.14) 



where <p is a smooth odd, increasing function on [— it, it] such that 

— , lor x > 

c/>(x) 



1 

-, tor x > — , 
2 2 

1 , vr 
, tor x < . 

2 2 



To be specific, we can choose 



<f)(x) = [rj* X[o,oo)](x) • X[-n,n](x) ~ 2' 

where \J is the characteristic function of the interval J and rj is a smooth nonnegative even mollifier 
such that 



/ ?7(x)dx = l,and \r](x)\ = 0, for all |x| > — . 
Jr 2 
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In order to keep H, identical with H] near ±7r, has to be near ±7r. We first note that there exists 
an odd function tp G (C%°(— n, n)) n such that (?/>,u ad ) = 1 since (C^°(— 7r, 7r)) n is dense in (L 2 (— 7r, 7r)) n 
and (u^,u ad ) = 1. We then define 

H 2 = Cj ijj(x - Oj), (2.15) 

where 

Cj = -(H), u ad (x - Oj)) - (Wj, u ad (x - - u ad (x)>. (2.16) 
Noting that Oj and Wj are small, this concludes the construction of Hj. 

Defining = {v G A c h | (vj,u ad ) = 0, for all j £ Z}, where X c h is defined in (2.5), we summarize 
the "smooth phase decomposition" procedure, denoted as <5phd 5 in the following lemma. 

Lemma 2.6 The "smooth phase decomposition" operator ^hd> as constructed above, is a smooth local 
diffeomorphism. More precisely, there are two neighborhoods of zero % G X c ^, "f G £ l x such that 
the nonlinear transformation 

^,hd: r — > ^ 

V = (0,W(x)) .— > {W^+H^ + u^x-^-u^x)}^ 
is invertible with ^ p hd an d =^hd smoo ^ rL - Its derivative at the origin is 



(2.17) 



Y=(0,W) ' — ► W + E*£, 
where E is defined in (2.18) below. 

Proof. We claim that 

(i) ^phd(O) = 0; 

(ii) 5p hd is C°°; 

(iii) ^h d (0), denoted as Jz? p hd, is an invertible bounded linear operator. 

Property (i) is straightforward. As for (ii), ^hd is smooth with respect to W due to the fact that ^hd 
is linear in for fixed 0. On the other hand, the smoothness of =5phd with respect to can be readily 
reduced to the smoothness of the mapping h-> {u*(x — Oj) — u*(x)}j e z- A direct calculation shows 
that, for given m G Z + , the mth-derivative mapping at is rj \- > {^ui (0j — x)r]j}j & z- We now only 
have to show that (iii) is true. In fact, the linear part of {H,- + u*(x — Oj) — u*(x)}j e z with respect to 
(0, W) around (0, 0) is E * = {Y,kez Ej~k0k}jez, where E = {Ej} jeZ with 

'^(x)-(| + i0(s)K(x), j = -l, 

-i(^(x) + ut(x)), J = 0, 

i^(x) - (i - i</)(x)K(x), j = l, 

0, others. 

Then we have the linear phase decomposition operator 

-£phd : x X^ — > X ch 

Y=(0,w) i — ► w + e*#. 



10 



Moreover, through direct calculation, it is not hard to obtain the bounded inverse of =S? p hd 

v i — ► (Fv,v-E*Fv), 



where 

^ch 



F : X, - ^ 



_^ r / U (2-19) 

By (i), (ii) and the inverse function theorem, the conclusion of the lemma follows. ■ 

Remark 2.7 The above lemma still holds when replacing X c h with ^ c h(XnH 2 ) and the proof is similar. 

In the new coordinates, the system contains lengthy expressions. We therefore introduce some simpli- 
fying notation first. 

6+ : C z — > C z 

X = {Xj} j( zZ ' > {Xj+l - Xj} j( zz- 



(2.20) 



5- : C z — ► C z 

X I > {Xj - Xj-l} j& . 

V : (C([-vr,7r],R n )) z — > R z 

X — > {(ViC-TT),^^))}^. 

Now, sorting out the linear terms, our lattice system is 

with boundary-matching and phase-decomposition conditions (2.11), (2.12), where N e / W represent the 
nonlinear terms of the system and 

A " = = ( id _ e , f) a - ( e * id ) = A ch - e. 5+ r) • (222 » 

where ^4 c h is the linear operator acting on the chopped variables; see (2.7). 

Remark 2.8 (i) Due to the fact that =^h/nf are isomorphisms, A n f = J^ p ^ d ^^ ^^h-^phd shares 
many properties with A. For example, A n { is sectorial in i 1 x X^ since A is sectorial in X. Here 
we use the definition of a sectorial operator from [12] which does not require the operator to have 
a dense domain. 

(ii) We relegate the detailed estimates of the nonlinear terms to Lemma 6.1 in the appendix since 
expressions are lengthy. We have, roughly, 

f |N e | ~ \(5+8) 2 \ + |0 3 ||5+0| + + |W|)(|W| + \5+d xx W\) + |W 2 | 

\ |n w | ~ \e\\6+e\ + + |w|)(|w| + \6 + d xx w\) + |w 2 |. 
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(Hi) Since the branch of continuous spectrum connected to A = may intersect the branches of con- 
tinuous spectrum in Re A < 0, it is in general not clear how to globally separate neutral from 
stable modes even linearly. Phase decompositions have been achieved globally in the case of weak 
pulse interaction, that is, in the regime where u*(a;) is close to a homoclinic orbit in the ordinary 
differential system Du xx + f(u) = ; see [18] for a linear analysis and [24] for a nonlinear 
reduction. 

3 Linear Four ier-B loch estimates 

In Section 3 and Section 4, we derive linear diffusive decay in our linear normal form 

V = A nf V. 

To illustrate the idea, we again use the linear heat equation 

u t (t, x) = Au(t, x). 

In order to obtain the diffusive decay on e Ai , we apply the Fourier transform and obtain the "diago- 
nalized" equation 

Ut(t,k) = —k 2 u(t,k). 

Then we have that \u(t, k)\ = e~ fe2 *|u(0, fc)|, for all t > and k G R, which, combined with Young's 
inequality, will give us diffusive decay for the scalar heat equation. 

In light of this procedure, we exploit Fourier transforms and the Bloch wave decomposition of A to 
construct an isomorphism diagram, from which we obtain a direct integral representation of A n { , that 
is, A n { = A n t (a)dcr. Unlike the explicit expression of e~ fc2 ', the estimates on e" 4 "*^* are more 

intricate and their derivation will occupy most of this section. 

To show the conjugacy between the linear normal form and its counterpart in a Fourier-Bloch space, 
we build a commutative isomorphism diagram involving the underlying spaces for these two operators, 
the linear operator A and its Bloch wave decomposition. To this end, we recall the definitions of the 
linearized operator A in (1.5), the chopping operator in (2-6), and the linear phase decomposition 
operator J2? p hd in (2.17) from above. We now consider these operators on L 2 /£ 2 -based spaces, that is, 
with new notation, 

A: (H 2 (R)) n — > (L 2 (R)) n , 

e 2 (Z,(L 2 (T 2lT )r) — ► (L 2 (m)) n , (3.1) 

i^ hd : e 2 x£ 2 ± (Z,(L 2 (T 27T )) n ) — > £ 2 (Z,(L 2 (T 2 ^)D, 
where T a = R/aZ is the one-dimensional torus of length a and 

£l(Z, (L 2 (T 2 ^)n = {w G £ 2 (Z, (L 2 (T 2n )) n ) | ( Wj ,u ad ) = 0, for all j G Z}. 
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(3.2) 



We write u = J" K u(x)e lkx dx and introduce several Fourier transform variants as follows: 

JF: I 2 — >■ L 2 (T!) 

i = {Qj}jez ' — ► J2jez d 3 e ~ l27rja > 

& n : (L 2 (HW)) n — ► (^) n 

u(x) i — ► u = {/^ r 7r u(j;)e- ifa dx} teZ , 

J^ ch : £ 2 (Z,(L 2 (T 2w )r) -> L 2 (T l5 (^n 

u(x) = {Uj^a u(a) = {E j6Z u J (x)e- i ^)(^'+-)}, 6Z , 

#- nf: £ 2 x^(Z,(L 2 (T 27r )r) — ► L 2 (Tx) xL 2 (T!,(£ 2 f) 

(0,u) T — ► (^(^,^ch(u)) T , 

where 

Ll(Ti, (£ 2 ) n ) = {w e L 2 (T!, (i 2 )") | ((w(a),^ n (e-^u ad ))) = 0, for all a G Ti}, 
We then have a commutative diagram of isomorphisms as follows, 

(L 2 (R)) n ^ £ 2 (Z,(L 2 (T 27r ))«) ^ 2 x^(Z,(L 2 (T 27r )n 

l^ 1 l^ch l^nf (3-3) 

Z/ 2 (Ti, (L 2 (T 27r )) n ) <^ L 2 (T l5 (£ 2 r) L 2 (Ti) x L]_(Ti, (i 2 )™), 

where is the inverse of the direct integral defined in (6.2), Section 6.2 and 

L 2 (Ti, (£ 2 ) n ) — > Z 2 (Ti, (I> 2 (T 27r )) n ) 

u(a) = { Ui (a)} jez .— > u(a) = (27r)^- 1 u((j), 

i£ hd : i 2 (Ti) x Zrj_(Ti, (£ 2 ) n ) — > L 2 (T!,(£ 2 r) 

(0(a), w(a)) i—). 0(<t)E(<t) + w(tr). 

Here we have 

1(a) = ^ch(E), (3.4) 
with E defined in (2.18). The inverse of Jz? p hd> which will be used later, has the expression 

J^J: L 2 (T 1 ,(£ 2 ) n ) — ► ^ ^(TOxLi^,^)™) 

w(cr) i — > (F(w(cr)),w(cr) - E(cj)F(w(cj))), 

where 

F: L 2 (Ti,^ 2 )") — > L 2 (Ti) 

w(a) — > -((w(a),^ n (e- ira u ad ))). 

We now use tildes for operators in physical space and hats for their conjugates in Fourier space. The 
index "ch" refers to the chopped operators, the index "phd" refers to the smooth phase decomposition 
operators, and the index "nf" refers to the normal form operators. We then define 

^ch := 1 ^4.=^ch) ^4ch := -i^chj (3-5) 
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where, according to the Bloch wave decomposition from Theorem 2 in the appendix, we have 

,1 

&T X A8 = A = \ B{a)da, (3.7) 



with B(a) defined in (1.7). Therefore, by the commutative diagram (3.3) and the equivalence relations 
in (3.6), (3.7), we find the conjugacy 

Anf = ^"f^nf^nf. (3.8) 

Just as we pointed out at the beginning of this section, based on this conjugacy, in order to obtain 
estimates on e Anft , we only need to derive estimates on e Anit . To this end, we first derive an explicit 
direct integral expression of A n f. From the equivalence relations in (3. 5), (3. 7), it is straightforward to 
see that 

A ch = [ 2 A ch (a)da, with A ch (a) := & n B(a)&~\ for all a G [-1/2, 1/2]. (3.9) 



Moreover, for any given (0(a), w(cr)) G L 2 (Ti) x L\(Ti, (£ 2 ) n ) and fixed a G [— \, \], by definition, we 
have, 



M w M = ^^A^ m L ) ) w 



' R(a) \/fl( ff )\ 

l ch (a)E(a) ^(cr) - |(<7)i2(<7) i lw(a)i 



(3.10) 



where 



We now conclude that 



F(a) : (l 2 ) n — )■ C 

w .— > -((w,^ n (e- ira u ad ))), 

i2(<r) : (^) n — )• C 

w i sillZ£(^(_l)V, J D< d (7r)). 

i . i 



(3.11) 



A ni = j 2 i A nl (a)da = f " ^ phd (<j)- l A ch (a)J? vhd (a)&o, (3.12) 



2 2 



where 

J? m (a): Cx(f 2 )«(a) — ► (* 2 )" 

(0,w) i — ► 6»E(a) + w. 

Here (^ 2 ) n (a) = {w G (f) n | ((w, ^„(e- ira u ad ))) = 0}. We also recall that A ch (a) is defined in (3.9) 
and E(cr) defined in (3.4). 



14 



Remark 3.1 We note that for any u E L 2 (T l ,(£ 2 ) n ) and v E 2t(A ch ), 

(Fu)(a) = F(a)u(a), (&6+T&£v)(*) = R(a)v(a), for a.e. a E [~, i]. 
In addition, for any (9, w) E £ 2 (Ti) x L 2 L (T 1 , (£ 2 ) n ), 



We now consider the family of linear systems, 

(^j = A ni (a) , for all a € [~ ±]. (3.13) 

While we obtained these operators based on L 2 /£ 2 spaces, we can also consider them on L q / £ q -based 
spaces. To be more precise, we first define a family of projections 

P q (°): Y q Y q 

w ^ w-^((w,^ n (e-^u ad )))^ n (e-^u'J, V ' ; 

where 

((£ q ) n , forls£g<oo, 
y g = T ' " 3.15 

[(^g°) n , forg = oo. 

Here we have = {x E | limini^oo |x n | = 0} with the supremum norm. For any q E [l,oo], the 
projection P q (a) is well-defined. In fact, ^ r n (e~ lax u a ^), ^ n (e~ iax u'+) E Y\ since u a d(±7r) = u^(±7r) = 0. 
We now denote Y q ^(a) = Rg P q (o~), and, in the following lemma, define A n f(a) on L q /£ 9 -based space. 

Lemma 3.2 For q E [1, oo] and a E [ — |, qL 

if phd (a): Cx7, iS ( ff ) — ► Y 9 

(0,w) i— »• #E(ct) + w, 

is uniformly bounded and invertible with its inverse 

^phd(^)- 1 : Y q — > Cx^M 

v i — >■ (-F(o-)v, v -^(cr)F(cr)v). 

Moreover, 

A; f ( ff ):Cx(y 9iS ( ff )nf ? (i>))) ->• Cx?,» 

is well-defined and sectorial. Here ^ q {A^{a)) = {w £ F 9 | {(1 + m- 2 )w m } m6 z E F g } is the domain of 
Ach(a) in Y q . 

Proof. The assertions for Jz? p hd(o") are straightforward. In order to show that A n f(a) is well-defined, 
we recall the definition of A n f(a) in (3.10), which indicates that we only need to show 

A ch (a)B(a) E Y q>s (a), Rg (A ch (a) - B(a)R(a)) C Y q>s {a). 
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We claim that A ch (a)E(a) G Rg(A ch (a) - E(a)R(a)). In fact, recall the definition of R{a) in (3.11) 
and define E{a,x) := {Y^ j E j (x)e- i2irja )e- iax G (C°°) n (^), we have 

R(a)E(a) = 2isin7Rj(E(cr, vr), £>u ad (vr)) = 0, 

which means that A ch (a)E(a) = {A ch (a) - E(a)R(a))(E{a)) G Rg{A ch (a) -E(a)R(a)). 

We now only have to show Rg (A^a) — E(a)R(a)) C Y gs (cr). Actually, for any w G & q {A c \ l {a)) with 
finitely many nonzero components, we have 

«^ch(<T)w -E( C r)i?((T)w,^ n (e- i ^u ad ))) =((l c h(^)w,^ n (e- ira u ad ))) + 27r J R(<r)w 

=2^(^(e iCT ^- 1 w), u ad ) + 27r(e iCT ^- 1 w, Du^O*))!-* 
=2vr(e ira ^- 1 w, J B*(0)u ad ) =0, 

and {w G ^ ? (A c h(cr))|w has finite many nonzero elements} is dense in ^ q {A c \ l {a)) under the graph 
norm of A c ^{a). Therefore, A n {(a) is well-defined. 

Next, A n f(a) is sectorial, due to the facts that A n ^{a) = ^f p hd(o") _1 ^ch( , 7)-^phd( ") an d ^ch(^) is 
sectorial (for details, see Section 6.3 in the appendix). ■ 

Now we are ready to obtain the estimates for the time evolution of system (3.13), for any given 
a G [— 5]- 0m discussion is split into the case a close to and the case a away from 0. 

For the case a ~ 0, the derivation of the estimate relies on a diagonalized normal form, that is, a 
complete separation of the netural and stable phase. First, we notice that spec(^4 n f (a)) = spec( J 4 c h(o")) 
is independent of the choice of q G [l,oo] and a G [—3,3], which we will prove in Proposition 6.4. 
Moreover, for a sufficiently small, there is a unique continuation of the eigenvalue 0, denoted as X(a). 
The set Ai := {X(a)} is a spectral set; see Section 6.4, 6.5 for detailed treatment. Hence, let 

Y q Y q 

be the spectral projection associated with A2 := spec( J 4 c h(o"))\{A(cr)}. Here e(a) (respectively, e*(cr)) 
is the eigenvector of the Bloch wave operator B(a) (respectively, the adjoint operator B*(a)) according 
to A(cr) with 

e(0) = ul, e*(0) = u ad , (e(a), e»> = 1. (3.17) 
We refer to Section 6.4 in the appendix for more details on e(tr) and e*(<r). We now denote 

y q ,c( a ) = span{e(cr)}, Y q , s (a) = Rg P q (cr), 
Ach(cT)\ YqA a) = Ma), A ch (a)\ Yq>s = A s (a). 
We then introduce the following diagonalized operator 



(3.18) 




It is not hard to conclude that for a sufficiently small, 

A dg (a) : C x (Y q , s (a) n ® q (A ch (a))) -> C x Y q , s (a) 
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is a well-defined operator. 

The key step here is to find an invertible bounded linear transformation 

% g (a) = ( T ^\ l 0l{ ° ] ] : C x Y g » -+ C x Y q>s (a) (3.20) 

such that ^ g (c")^4 n f (a) = A dg (a) ^g(c) . We note that the choice of =^d g (cr) is not unique since there 
are nontrivial invertible operators that commute with A dg (a). 

Lemma 3.3 For a sufficiently small (that is, \o~\ ^ jq) and q G [l,oo], 

^ ~ \p q {am°) W\y qA J ( } 

satisfies the relation =^d g (cr)^4 n f (a) = Adg(a) ^d g (o~) . Here we have that fi(a) = — ^:((E(o"),«^" n (e*(o")))) 
and 

S(a): Y q — > C 

w ' ^ -^((w,^„(e*(a))». 

Moreover, we have 

^dg 1 = (^00 — ^oi^i/^io) 1 f /fi-i/fi ,ffi ffi rp-lrr yr-1 \ /f^-l/f? m /f?-l ) ' (3.22) 

\ _J ii J io Uoo - ^ oi-t ii -tioj-tn + 1 n J io J oi J n y 

in which we suppress a -dependence for simplicity. 
Proof. We recall from (3.10) that 

^=(id-ItW))^ w ( £w ,d )- 

Therefore, in order to find a ^ g as required, we only need to find an invertible bounded linear operator 

^nt(a) = ( :Y q ^Cx Y q , s (a) 

such that 

^nt{o-)A ch (a) = A dg {a)^ nt {a), 

which is equivalent to 

\^(a)A ch (a) - A ch (a)^(a) = 0. 

While the choice of =^1/2(0") satisfying the above equation is apparently not unique, we choose that 
Ji(ct) = 5(a) and = Pq{o~)- As a result, we have 



£ W _W fE(a) idl- )- M( 5 
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To show that ^d g ((r) 1 in (3.22) is correct, we only need to verify that 



which is clearly true. 



Based on this lemma, we now derive the estimate for e A ^(°') t when a is close to zero. We first introduce 
new notation M(t,a) := e Anit and J%(t) := e Anit with 

M<* ^ (M a(t,a) M x(t,a)\ ~ ~ 

/ \ (o.2oJ 

^ (<)= Ur"(t) 1h(!)) : ^^( Z '( i2 M) ~> £ 2 x£i(Z,(^(T 27r ))-). 

To make sense of the derivatives and Taylor expansions with respect to a of entries in ^ g (o~), we 
extend Toi(<r) and Tn(o~) continuously as operators on Y q , that is, 

f 01 (a) = S(a)P q (a), f n (a) = P q (a)P q (a). (3.24) 

The same argument applies to operators ^g 1 ^) and M(t,o~). 

Lemma 3.4 For a sufficiently small (that is, \o~\ ^ 70J and q £ [l,oo], there exist positive constants 
C(q) and d such that, for all t ^ 0, 

{\\\M 10 (t,a)\\\ c ^ Yq \\\Mn(t,a)\\\ Yq ) - ^ & J ' ( > 

Moreover, we have a higher regularity result for M\\{t, a), that is, for any given a G [— 70, 70] , q £ [1, 00] 
and a > 0, there exists C(q,a) > such that 

\\\Mn(t,a)\\\ Yq ^ Yqa < C(q,a)[(l+t- a )e-^ + 3-^6"**]. 

Proof. The idea is to evaluate M(i, cr) = e^» f ( ff)i based on e^W* = ^ g (o-) _1 e^ d s( CT )* J^ g (cj). We first 
state the following estimate from Proposition 6.7: For all q E [l,oo] and a £ [—70,70]) there exists a 
constant C(g) > such that 

|e A ^| C(q)e-$° H , \\\e^% < C^e"**. 

To obtain estimates on ^ g and its inverse, we start by computing the Taylor expansions of entries in 
^ g (o"). A straightforward calculation using (3.17), (3.4) and (2.18) shows that 

e(o-) = < + icrei + 0(a 2 ), e*(a) = u ad + iae\ + 0(a 2 ), 

e -iax = 1 _ i(JX + o(a 2 ), E(o-) = - 27rio-0< + io-xu'J + O(o- 2 ), 

where ei (respectively, e^) is even and nonzero due to the fact that B(0)ei = — 2Du" (respectively, 
B*(0)e\ = —2Du' ad ). Then, plugging these expansions into 5d g , and using (3.24), we obtain 

~ (tT s _ ( \ _ / 1 -2™*\ ( 0(a 2 ) 0{a 2 )\ 

^ [(T) U(a)E(a) P q (a)P q (a) \-2ri** n (*) P q (0) + \o(a 2 ) O(a) I' ^ 
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where = — + </>u* — and 



9 

w 



In 



^((w^nCxUad + e;))). 



Therefore, for any cr E [— 70,70]) <7 G [1>°°], there exist positive constants C(q) and d such that, for all 
(6, w) £Cx Y g , we have the following estimate. 



\M 00 (t,a)9\ \M 01 (t,a)w\ \ < / 1 |cr|\ /e^ 
\M 10 (t,a)6\\ Yq ||Mn(t,(7)w||y J " W \M 1 







e 2 



1 H 
kl 1 



< C(g)e 



1 I erl \ _d 2, / lar I cr I \ 
M kP 6 2 + J| 1 6 2 



|0| 



Using (3.20), (3.22) and (3.23), we now expand M%i(t,a) and obtain 

UMiifoaJWIIya ^C(||fi 1 (a)- 1 f 10 (a)||Y 9 c,|e A ( CT )*foi( ( 7)W| + \\f n (a)- 1 e^ t T 11 (a)W\\ Yr + 
llff/^fioC^llvalfoi^fi^^-^^^fnCa)^!) 



<C(g, Q)[|a| 2 e-f CT2t + + lje"^* + |cr| 2 e 

. 71 



|W| 



^C(q, a)[(t~ a + l)e-^ + ^-e^ 2 



l + t 



at ]\m\Y a 



where in the second inequality we used (3.26), and Proposition 6.7. ■ 

Remark 3.5 We point out that, in the above lemma, the estimate for Mio(cr) can not be improved, 
since #" n ($) 7^ 0. In fact, due to the fact that <£(x) G (C°°(T2 7r )) n and is a nonzero even function, 
we have 

B (0)Z> = ^[-2Z)< - 5(0)6! + 2ttS(0)(K)] = B(0)(<f><) + 0. 

On £/ie other hand, the estimate for Mq\ (cr) can be improved given suitable additional assumptions. For 
example, if we assume that u' d (±7r) = 0, then a;u ac j + e\ is zero, which leads to a better estimate. 

For the case cr away from 0, we have the following result. 

Lemma 3.6 For a away from zero (i.e., for 70 ^ \o~\ ^ \) and q 6 [l,oo], there exist constants 
C(c/)> 72 > such that 

\\\M(t,a)\\\ CxYq ^C(q)e^ t (3.27) 

Moreover, we also have a higher regularity estimate for Mn(t,a), that is, for any given 70 ^ \o~\ ^ \, 
q G [l,oo] and a > 0, there exists C(q,a) > such that 



\M n (t, o-)\\\ Yq ^ Yqa < C(q, a)(l + t~ a )e 



Proof. Recall that e A ^) 1 = ^^(a)' 1 A^' : -£?phd(<r) • The inequality (3.27) is true due to the 
uniform boundedness of ^(o~) in Lemma 3.2 and the fact that |||e j4ch '- cr )* |||y ^ C(c/)e -72 *, for a away 
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from 0, in Proposition 6.7. Moreover, by the expressions of J^hdC ") an d its inverse in Lemma 3.2, we 
have M u (t,a) = (id - S(cr)F(cr))e j4ch(<j)< . Applying Proposition 6.7, we conclude that 

\\\Mn(t,a)\\\ Yq _ Yqa = HI (id -l(a)F(a))e^W*|||^ Q < C(g, a)(l + r")^*. 



Lemma 3.4 and 3.6 give the following proposition. 

Proposition 3.7 (Fourier-Bloch estimates) For any a G [— |, |], G [l 5 oo], i/iere exzsi constants 
C{q), c > suc/j i/iai v4 n f(cr) is sectorial and 



( \M 00 (t,a)\ |||M 01 (t,cT)||| y ^ c \ / 1 



\\\\M w (t,a)\\\ c ^ Yq \\\M u (t,a)\\\ Yq 



< C{q) i e ~ CU > /° r a// * > °- ( 3 - 28 ) 



/f+1 *+l 



Moreover, we have a higher regularity estimate on Mn(t,a), that is, for any a G [— h, h], q G [l,oo] 
and a > 0, there exist constants C(q,a), 7 > such that 

\\\Mii(t,a)\\\ Yq _>y* ^C(q,a)(^(l + t~ a )e-^ + -^e- ! l a2t ^, for all t > 0. (3.29) 

We also need the Fourier-Bloch estimates for the derivative d a M(t, a) in the following lemma. 

Proposition 3.8 (Fourier-Bloch estimates for derivatives) For any a G [— hih]> Q. £ [1>°°] an d 
(3 G (2(1 — -), 1), there exist positive constants C(q, (3) and c such that, for all t ^ 0, 



\(d a M) 00 (t,a)\ |||(^M)oi(*,a)||| n _ c 
|(0 ff M) lo (t,<T)||| c _ yo |||(5, 



^ff^) * ( ' ^+^ P >~^- (3-30) 

^M) n (t,a)\\\ Yq J j±r t J 



Moreover, we have a higher regularity estimate on (d a M)n(t,a), that is, for a G [—§,5], q G [l,oo], 
P G (|(1 — -), 1) and a G (0, 1), there exist C(q, a, /3) > and 7 > such that 

\\\(d a M)n(t,v)\\\ Yr + Yqa ^C{q,a,P) ~ ^ aH + (^" tt + t 1 ^)e~^ , for all t > 0. 

Proof. On the one hand, we take the partial derivative of the following system with respect to a 

\mt,v) \W(0,a)l 



and obtain 



A.) = M(t '" ) UxM + (8 » M(^ "' ), W(0,a) 



On the other hand, we have that 
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Taking the partial derivative with respect to a, the equation becomes 

Uw)(«,<J = Uw (t , J +A "< (<r) 



for which the variation of constant formula gives 

U„W(t,<7)( v ' ' l8„W(0,<r)J y„ V \.W(0,<rW 



Therefore, one has 



d a M(t, a) = [ M(t - s, a)A' nf (a)M(s, a)ds 
Jo 

J4M(a)-W-)(*-*)j4 M (a)^ (3.31) 
^ P hd(<T)- 1 e^ h(f7)( *- s) ^( C T)e^M^ phd ( C 7)d S , 



where 

^(<r) = A' ch (a) + A ch (a)E(a)F(a) - E (a)F(a)A ch (a). (3.32) 
We recall that E_(a) is defined in (3.4), F(a) in (3.11) and A ch (a) in (3.9). 
For \a\ ^ 70, by Lemma 3.3 and the above equation (3.31), we have 

d a M(t,a)= f ^ g (a)~ 1 JF(a,t,s)^ g (a)ds, (3.33) 
Jo 



where 



JK(a,t,s) = e^WC*-) f^j (-^„e(a), id) e**M- 

_ / -e A (' J )*S((j)^K((j)^ n e( f 7) e A M(*- s )5(a)^(a)e^ s 
" l v _e A M s e ls M(*- s )p g (a)^(a)^ n e(CT) e^^P^a)^^)© 1 ^ 

/ iVoo(cr,t,s) iVoi(cr,t, s) \ 



(3.34) 



' \N w (a,t,s) N 11 (a,t,s)J ' 

We now evaluate the entries of JV with expansions combining (3.32) and (3.34). First, recall the 
definitions of Ad g (cr) in (3.9), S(a) in Lemma 3.3, P q (cr) in (3.16), in (3.2), and e(<r) in (3.17). 

For -/Voo, note that it is smooth with respect to a and 

N 00 (0,t,s) = -5(0) (A' ch (Q)+A ck (Q)M(0)F(0) -E(0)F(0)A ch (0)) J^e(0). 



We claim that N 00 (0,t,s) = 0. In fact, since ^4^(0)^(6(0)) and A ch (0)E_(0) are orthogonal to 
^ n (u ad ) in Y 2 , 5(0)^(0)^,(6(0)) = and 5(0)! ch (0)!'(0) = 0. Moreover, F(a)A ch (a) = R(a), 
which is defined in (3.11) with R(0) = 0. Therefore, there exists a positive constant C such that 

|iV 00 | ^ C\a\e~^ 2t . 
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For Niq, due to Proposition 6.7 and the fact that Niq is smooth in a with iVio(0, t, s) 7^ 0, there exists 
a positive constant C such that 

lll#iolllc->y,, W < Ce-^e-l^ < Ce-^'e-^C*-). 
For Nqi, we have, for any q G [1, 00] and /3 > ^(1 — |), 

lll^i||| y „. (<rH c (111^(^(^111^^+ 

|5'( t T)^ ch ( C T)^ / (<T)||||e^^> s ||| n>s(tT) + |||i ? (cr)A ch ( f T)e^^) s ||| nis{CT) ^ {C ) 
^Ce-s^C*--) f llle^^lll 1 + klllle^^lllv , a + M\\\e A ^> 



^C(/3)e"2 CT ^ \\\e A ^ )s \\\ 1 + |||e As ^ s ||| y M + |<r| |||e A ^ s ||| y ,, Y ? , 

where the last inequality results from the fact that, for any g G [1, 00] and /3 > ^(1 — -), we have a 
continuous imbedding 

Now, using Proposition 6.7 and 6.5, we can further conclude that 

lll^>l|||y ais(CTHC **C{q,P){8-* + 1 + \a\s-P)e-^ 2 ^e-^ 
^C{q,^)e-^ 2t (s-^ + \a\ + H^e"^. 

For iVn, we have, for any q G [1, 00] and /3 > |(1 — ~), 

III#ii|IIy..m ^^lll els(,7)(t " s) llly os M (lll^)e^ 



,( ff ) ^iif V 111 chV ; mY q;S (<r)^Y q 

+ \\\^ )s \\\ YqAf7) + ll|F(a)^c h (^)e ls(CT)s ||| yg , s(CTHC 



Therefore, combining (3.33), (3.26), and the above estimates for entries, we conclude that, for |<r| G 
[—70,70], q G [l,oo] and j3 G (^(1 — |), 1), there exist positive constants C(q,(3) and c\ ^ | such that 

/ I (a CT M)oo (t, a)| HI (daM)oi (t,a)\\\ Y ^ c \ 
{\\\(d a M) 10 (t,a)\\\ c ^ Yq IKd^u^a)^ ) 

1 H\ f* ( Woo(v,t,s)\ lll^oi(<7,M)|||y 9iS ( CTH c 







W,:, 







1+t 







Hi 





(3.35) 



1 1 



Here the inequality (*) relies on the fact that for any /3 G (0, 1), there exists a positive constant C(/3) 
such that 

-* - r* „ -n.. _ t 1 -^ 



^ e^'da < C(/3)--^, jf S -V>d S < COS) 



+ 
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On the other hand, for 70 ^ \a\ ^ |, q G [1, 00] and (3 G (^(1 — |), 1), by the expression (3.31) and 
Proposition 6.7, there exist positive constants C(q,f3) and C2 such that 



\\\d a M(t,a)\\\ CxYq jf |||e A -W( i - s )^( CT )e^M s ||| g d S 

<C(g, /3)e" 72 ' / (a"3 + 1 + |a|s- /3 )ds (3 - 36) 
Jo 



"21 2]' ? G L J 



By (3.35) and (3.36),we now conclude that, for any cr G [— A, A], q G [l,oo] and j3 G (|(1 — ^), 1), there 



exists positive constant C(q,j3) and c such that 

/ KMOooft*)! lll(^M) 01 (i, CT )||| y ^ c \ / 1 ^\ 1 2t 

^||iaM) 10 (t,a)||| c ^ |||(ftr^)ll(t,0-)|||y, ; ^ \7fe A/ 

We now consider (<9 CT M)ii(i, a). For a G [— 7o,7o]> we P m g =^dg( cr ) from (3.20), ^g(o-) -1 from (3.22), 
and N(a,t,s) from the last equality in (3.34) into (3.33). We then obtain 



(d a M)ii(t,cr) 



ds 



J -ff/fxo (iVooToi + iVoifn) + ff/ (iViof i + iVufn) ds • (1 + 0(a)) 
More precisely, for q G [1, 00] and a G (0, 1), there exists C such that 

\\(d a M)u(t,a)W\\ Yq ° J [ife^iolly- (|iV 00 ||T 01 W| + |iV if n W| 

+ ||f n 1 iV 10 ||y (r |fo 1 W| + llff/iVnfnWHy. 

Zc jf [||fi ||ya (|iV o||?oiW| + |^oiT\iW| 
+||^io||K»|f iW| + ||iViiTiiW||y«j ds. 
Here the inequality (**) relies on the fact that 

T n (<7): Y q * — >• F« 

v ' ^ v-(e(a),e- i ^u ad )- 1 ((v,^ n (e- i ^u ad )))^ n e(a) 

is a uniformly bounded operator for (/ G [l,oo] and a G (0,1). Using the explicit expressions of the 
entries of ^ g (cr) in (3.21) and the estimates on the entries of jY as shown above, we derive the following 
estimates, 



[ ||fio||v-|JVbo||r iW|d« < C(g,a)|a| 3 te-^ 2 *||W||y 9 , 
J 

J* \\T 10 \\ Y ~ |iVoif u W|d S < C(q,a,(3)\a\ + H (1 ^^_J e-^ 2 *||W|| n 



,t 1 



/ ||iV 10 ||y«|T 01 W|d S < C( (7 ,a)^L=e-f CT2 *||W|| y9 , 
jf HiVnfnWHy.ds ^ C(q, a, 0) ft + t 1 ^ + J^{t 



i 
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We now conclude that, for a G [— 7o,7oL 1 £ ft £ — f)>l) an d a ^ (0, 1), there exists 

C(q,a,(3) > such that 

i 1 _ /3 

\\\(d„M) u (t,a)\\\ Yq ^ Yqa <C( q ,a,P)[ t2 + + t 1 ^)^}. 

For 70 ^ |cr| ^ |, g G [1, oo], /3 G (|(1 — |),1) and a G (0,1), there exists C(q,a,(3) > such that 
||P CT Af) n (i,a)||| Yg ^=|ll(id -l(a)F(a))e^('K^)^( ff )e^W*||| y ^ a 

^(7(g,a,/3)e- 72 * / (t - s)- a s"5 + s "3 + 1 + \a\s^^ds 
Jo 

^C(g,a,/3)(t5- Q + i 1 ^) e -^. 

Altogether, for a G [-|, |], q G [l,oo], /3 G (|(1 - |), 1) and a G (0, 1), there exist C(q,a,/3) > and 
7 > such that 

- 1-/3 

HI (d a M) n (t, a)\\\ Yq ^ Y » ^C(q,a,P)[ t2 +* t e~^ 2f + (^" a + i 1 "^"^], for all t > 0. 



4 Linear estimates in physical space 

According to the outline at the beginning of Section 3, we are now ready to derive the linear estimates 
for e Anfi . To be more precise, we first show by Fubini's Theorem that 

[^j = M(t) * , 

where M{t) is the generalized "inverse Fourier transform" of M(t,a). We then employ an argument 
similar to, but more intricate than, Young's inequality for the case of the scalar heat equation, exploiting 
the linear Fourier-Bloch estimates in Proposition 3.7 and 3.8, to obtain the general L p -L q estimate on 
our linear normal form e An{t . 



To this end, we first note that A n { = A n f\p ly , x ± and thus we have, by (3.8), for any (9, W) G i 1 x X^, 

■*® (ffi) " c """ (l) " """" (w) - (|) • for ■» ' > °- 

Recall the notation ^(t) = e An{t , the definition of J^ n f from (3.2), and the definition of A n f, A n { from 
(3.6). In addition, by (3.12), we have, for any (6(a), W(<r)) G L 2 (T X ) x L^Ti, £ 2 ), 



c 
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To show that e Anft is a generalized convolution, we first define M(t, cr)'s "generalized inverse Fourier 

transform" M(t) := — 00 — 01 , with expressions as follows 
yMio Mn J 

i 



M 00 (t) :={M o(t,i)b 6 Z := { / " M 00 (t,a)e i2 ^da} j6Z , 
M 01 (t,y) : ={M i(t,y,j)},a := { £ £(M„iMt, a)e-^y^da} jeZ , 



"5 tez 
i 



M 10 (M) :={Mi (t,x,j)} ieZ := X)(MioM*, C r)e i ^ a; e i27r ^da} ieZ , 

1 

Mu(t,x,y) :={Mn(t,x,y,j)} ieZ := {— £ £ (Mn^f^je^^^e^+^e^^da} 



2 f!,r/e2 



We then have the following lemma. 

Lemma 4.1 For any (9,W) G i 1 x and a// 1 > 0, 



where 



MOO * £ ={J] Moo(*, j " fc)0 fc }i6Z, 

fcez 

Moi(t,y,j- k)W k (y)dy} j& , 

-7T 



MlO * £ ={J^ x, j - fc)0 fc } jeZ , 



(4.1) 



wl w IwJ \Mio*i K u *w v ; 



fcez 

Mn(t,x,y,i-fc)W fe (y)dy} i6Z . 

Proof. The proof is a straightforward application of Fubini's theorem. ■ 

We are now ready to obtain the general LP — L q linear estimates on ^K{t). We denote 

X q = (^(Z,^(T 27r )))", for any q G [l,oo], 
and prove the following proposition. 

Proposition 4.2 (general L p —L q estimates) For any 1 ^ q ^ p ^ oo and (9, W) £ x X^, £/iere 
exists a positive constant C such that, for all t > 0, 

/ii^oo(^ii^ ii^oi(t)wii*\ < / (l+tr^i-^ngi^ (i+t)-^-^\\w\\x g \ U3) 

\\\^io(t)e\\ Xp \\^n(t)W\\x p ) " ^i + t )-|(i-?)-i||^||, 9 r^i-F)(l+t)-i||W|| x J' 
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Proof. We illustrate the derivation of the estimates on ^#01 an d sketch the estimates on ^#00 an d 
^#10- Lastly, we show the estimates for 

We first notice that, for any W G and 1 ^ q,r ^ p ^ 00 satisfying 1 + - = | + ^ , there exists a 
positive constant C such that 



IKoi(t)W||^ ^ C||M, 



1 1 

9 



M J] sup |M i(t,y,j)| ||W|| Xg . 



(4.4) 



In fact, by Holder's inequality, we have 



oi(i)W||f p =||M 01 *W||^ = ^|^ / M Q1 (t,y,j-k)W k (y)dy\i 



^llMoiWIirilWlir/ £ f l*oi(*.y.J-*)riW fc (y)|«dy 



j,kez 



^\\K i(t)\\ P x - r 



sup^|M i(t,y,j)r ||W||^ 



(2vr) 1 "l||M 01 (t)|||J [ sup ^|M i(i,y,j)| ] ||W| 



-V„ 



Moreover, by (4.1), we have 



IMoiWIIx 



sup I / I Y,{MoMt,o)e-^y\do\ < — ^ 



C(oo) 



2 £e 



1 + t 



(4.5) 



Here we use the fact that any bounded linear functional on £q° can be viewed as a bounded linear 
functional on £°° with the same norm. We now estimate the X\ norm of {Moi(t,y,j)}j^%. By using 
Proposition 3.8, there exists C > 0, independent of the choice of y G [— 7r,7r], such that 



E l*oi(t, y, j)l =e( 1 + ~" (l + { -LJA_) 



\M 0l (t,y,j)\ 



7— 

Jk 1 + ^ 



-dx 



E( 1+ (J / )a )l^bi(t,y,j)l 



I 1 

2 



2 a=0 £GZ 



*** 1 / / 2 g 



1 |-f ( t | + t l-f)2 e -2«T a t 

TTt 



— d(T + - , 

_I l + t t J_l 



da 



(4.6) 



u + 1 



(i + t)l ^/TTt 



for all i > 0. 
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Here in the inequality (***) ; we applied Proposition 3.8 with q = oo and /3 = | (actually, any fixed 
P e Combining (4.4), (4.5), and (4.6), we have that, for all 1 ^ q ^ p ^ oo and € -^^j 

there exists a positive constant C such that 

C 



\\VT+t^ 01 {t)W\\ eP ^ _^^||W||jf„ for all t > 0. 

(1 +t) al 9 f J 

For ^#00) the steps are the same as above but easier. For ^#io, we point out two main differences to 
the above calculation. First, instead of (4.4), we use 



\\^Mm\\x p < C\\M 10 (t)fxJ 



J 7 " l^\M w (t,xJ)\ ) dx 



V 



Second, to estimate the Y\ norm of {Mio(t, x, j)}jez, we use Proposition 3.8 with q = 1 and f5 
^(actually, any fixed (3 £ (0, |]), instead of q = oo and = | . 

The last step of the proof consists of deriving the estimates for dt\\. We first have 



i_ i 

q p 



\\^u(t)W\\x p ^(27r)r sup sup \M u (t,x,y, ^p > \M u (t ,x,y,j)\ ||W|| 

\|x|,|y|<K|i6Z / \\x\,\v\<«5£, J 

On the one hand, we apply Proposition 3.7 with q = oo and a > ^ and have 

sup suplMufrx^OK sup / 2 | y2{M u ) ev (t,a)e iia+e)x e- [{a+v)y \da 
|x|,|y|<K|i6Z M.lvKM-M ^ z 



^C(a) / i |||M 11 (t,a)||| yoo ^d C 7 
C(a) 



t a (l + t)2~ a 

On the other hand, by applying Proposition 3.7 and 3.8 with q = oo, a E (i, 1) and /3 = |, there exists 
C(a), independent of choices of x,y £ [— it, it], such that 



_i i 

x~y\2 \ 2 / /■ | x-y\2\ 2 



J2\Mn(t,x,y,j)\ = J2 [l+ ij + y ) J 2 h + (j + ^ ) |A/ n (t,x,y,j)| 



lil>i b>i 



.1 i 



<C(a)** N ° |||Mn(t, ^IH^^da + £ HK^M)^, ^W^^da 



<C(a)t 



2 2 
1 
2 



1 



^ 2a (l + 



<C(a) , = — , for all t > 0. 
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Moreover, combining the above two estimates, we have that, for given a G (|, 1), there exists C(a) > 
such that 

C{a) 



sup ^jMnfoa.y.j)! < 

r . In Kit . _™ ' 



Therefore, for any 1 ^ q ^ p ^ oo, a G (^, 1) and G -X^,, there exists C(a) > such that 

ii^uwwiu, « (i+ ( ; ) ( °L Ha(1+ 1 t)1 _j iwiu,. 

Moreover, we can improve the above estimate for t close to zero. Note that for the Laplacian operator, 
we have the general L p -L q estimate for all t > 0. As a perturbation of the Laplacian operator, has 
the same estimate for sufficiently small t, which can be seen by using the variation of constant formula 
as follows. 

Il-^n (*)W|| Xp = || (id -E*F)e^*W|| Xp < C\\e A ^W\\xp = C||e A 'W|| L p, 

where W = {Wj(x)}j & z and W(2ttj + x) = MVj(x) for all j £ Z and x G [— 7r,7r]. We now let 
V(t, x) = e^W(i) and have 



V(t)=e %I 'W+ I e D9 ^ s) f'(u ir )V(s)ds. 
Jo 



from which we derive 

sup ^(?-|)||v(t)|| iP ^ ||W|| L 9 +CT 1_ ^ ( 9 _ p ) sup ^(?-i)||v 

0<t^T 0<i^T 

1— -(---) i 

Taking T sufficiently small such that CT 2( i p' ^ ^, we obtain 

||e A *W|| L p ^ f_ x ||W|| L ,, 

f2(q p> 

which implies that 

\\^u(t)W\\x p < rr l|W|| Xa , for all < t < T. 

Therefore, for any 1 ^ q ^ p ^ 00, there exists C > such that 

C 1 



LP- 



||ufn(t)W||x„ < lf i 1W1 ,,J |W||x, 



|) (1+t)' 



Remark 4.3 5?/ (4.1), (4. 2) and a similar argument as in Proposition it is not hard to conclude 
that, for any j G Z + , 1 ^ q ^ p ^ 00 and (9,~W_) G ^ x -AT^, f/iere exists a positive constant C such 
that, for all t > 0, 

/ ||<5^bo(t)£||* |l4^oi(t)W||^A 
l|l4^io(*)£||x p Il4^ii(*)w|uj 



^ c 



(l + t)~^-p +j) \\0\\ eq (l + t)-i { ~p +3+1) \\W\\x q \ 

(i + ty^-p +j+1) \\e\\ £q t-* ( *-p\i + t)-( 1+ i^w\\ x J ' 
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5 Maximal regularity and nonlinear stability 



In this section, we prove the main theorem-Theorem 1. To achieve this, we first introduce a Banach 
space that our argument will be based on. We then collect several maximal regularity results since the 
normal form system is quasilinear. Based on our normal form and the general L p — L q linear estimates, 
we can apply a fixed point argument to the variation of constant formula, thus obtaining the nonlinear 
stability result. 

We choose r G (4, +oo) and define 

z = {(0, w) e c((o,+oo),^ x (x ch n ^ ch \H 2 ))) | YOh < °o}, 

where 

\\(9,W)\\z =sup \\6(t)\\ e i +sup(l + ^ \\9(t)\U™ +sup(l + t)i\\d 2 ey 

+ sup(l + i)2 ||W||xi + sup(l + 4)|| WUxoc + + sup(l + t)^||<5+W||x a 

i>0 i>0 i>0 

1/r 



+ y™(i + ty\\5 + d xx w(t)\\ r x?) 



Here we have 5 2 := 5-5+, where 5± is defined in (2.20). 

Lemma 5.1 (maximal regularity) For any given T > and r £ (l,+oo), there exists a positive 
constant C such that the following holds. If (77, v) G L r ((0,T),£ 2 x X2) and if (9, w) satisfies 



'9_{tY 
w(t) 



j^(t-s)^^ds, te[0,T], 



then 



J xxl 



w(t)\\ r x 2 dt ^ C (|| 2 (*)ll<- + l|v(t)l|jc a ) r d*- 



Proof. The result just follows from the standard maximal regularity results on the Laplacian operator 
and the robustness of maximal regularity with respect to lower order perturbations. To see that, we 
first recall ^(t) = e nft , where A n f is defined in (2.22). By [13], the maximal regularity result holds 
when we just replace A n { by Aq, which is defined as 



A) = 

Viewing A n { as a perturbation of Aq, we have 

rt /„/„\\ rt 





Dd xx 



Then by the maximal regularity property of Aq, we obtain 

rT rT 



L 



^ \\d xx w{t)\\x^ ^C[ T (ll^nf - A ) \\ fixX2 + II \\ PxX2 \ ds. 
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We observe that, for any e > 0, there exists K(e) > such that 



||(Anf - Aq) I - J \\ £ 2 XX2 < e\\A I - I y xX2 + K(e)\\ I - I || £2xX2 . 



In addition, it is straightforward to see that 

The conclusion follows by combing the above three inequalities and taking e sufficiently small. 

We also prove a corollary which will be useful in the proof of nonlinear stability. 
Corollary 5.2 For given a£l and r G (l,oo), there exists a positive constant C such that, if 

r-t 



then 



poo poo 

1 (i+tr\\d xxYL (tw X2 dt^c (i+tr {{my + \\Y.(t)\\x 2 ) r dt. 

1 JO 



Proof. We first note that, for t G [n,n + 1), n G N\{0} 



w(t) ; vA-i L-i y v x(s) i 



r/(n — 1 + s) 



- + 1 - s) - ; ( ds. 

\ v(n - 1 + sW 



rt— n+l ft—n^ 
10 JO 

Applying Lemma 5.1 to the above expression, we obtain 

'n-1 

The conclusion follows from multiplying both sides with n a ~ (1 + t) a and summing over n G N\{0}. 



/ ||0xxw(t)||5r a d< < C / (||r?(t)||, 2 + ||v(t)|| X2 ) r dt. 

Jn J n—1 



Lemma 5.3 If (0q, W ) G i? 1 x (X^ n (ii~ 2 )), £/ie solution of the linear system 
belongs to Z and there exists a positive constant C\ > suc/i i/iai 



W(t) j " Z ^ °^ I Wq J l^x^n^ 1 ^)) 



^^IIUxt- lUxOcn^rVi^n- ( 5 - 1 ) 



30 



Proof. By Proposition 4.2, it is straightforward to see that 

\\^oo(t)e \\ el < CWioWtL, ||^bi(*)Wollp < (1+ ^ )1/2 l|Wollxi, 

c c 

||^bo(*)foll*» < / 1 + t u /a l|goll^» ll^oi(t)W ||^ < -_||w || Xl , 

c c 

||5 2 ^oo(^oll^ < (1 + t x 5/4 l|goll^, P 2 ^oi(i)W || £2 < ^J IWqIUi, 



IKIO^IIX! < (1 + t) i/2 ll^ll^' H^llWWolUx < ^llWollXi, 

ll^io^ollx^ < -^r t \\e y, Mn(t)w \\ Xoo < J^||w || Xoo . 



Moreover, we have 



\5+d xx W(t)\\x 2 ^ C\\6 + A ^(t) (~ u I y xXa 



< C(\\5 + A ni ^(t) I -J ||, 2xX2 + p + (A nf - A))^T(t) L^J ||^xjr a )- 

We need to show that the two terms on the right hand side of the above inequality decay sufficiently 
fast. On the one hand, we claim that 

\\s+A ni ^(t) (jA y xXa ^ ^^rll fel H^x^ 1 ^). fOT a11 * > o. 

Actually, for t G [0, 1], the above inequality is true since 5+ is bounded and 

\\A at ^(t) (jjj^J \\ PxX2 < Cp nf f^J ||,2 XX2 < C|| ^ || £2x ^ 1( ^ 2) . 

For f £ [ljOo], we first point out that, to show \\\A n f^(t) = A n fe Antt \\\ i 2 xX2 decays with rate i _1 
as t goes to oo, we only have to show that the supremum norm of its Fourier-Bloch counterpart 
A n f(a)M(t, a) decays with rate t~ l as t goes to oo, just as in the scalar heat equation case. This is true 
by applying the steps in Lemma 3.4 and Lemma 3.6 to A n f(a)M(t, a). Second, it is straightforward to 
see that the discrete derivative operator <5+ gives an extra t -1 / 2 decay, which concludes our justification. 
On the other hand, we have the explicit expression, using that <5 + and A n f, Aq commute, 



6+(A n{ - A )^f(t) 



110 



S+T \ 5 + 8{t) 



v w y \A ch B* f(u*)-E*$+r/ \6+W(t) / 
We apply Proposition 4.2 again and obtain 



||A ch E * (S+9(t))\\ P < C\\5 2 6(t)y < — ^-r (ll£oll^ + l|W IU 2 ) , 

(1 + *) 7 



4 



|f'K)(5+w(t))||x 2 < c\\s+w(t)\\x a < —^(U y + l|w IU- 2 ). 
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In addition, recalling that T is defined in (2.20), we conclude that, for any e > 0, there exists K(e) > 
such that 

\\6+T(6 + W(t))\\ Xa < C\\d x 6+W(t)\\ Xa < e\\5 + d xx W(t)\\x 2 +K(e)\\5 + W(t)\\x 2 . 
Therefore, by choosing e sufficiently small, we conclude that 

||* + M W(t)||x a < — ^-r (114,11/1 + l|W ||^ 1(H2) ) , 



which shows that 

l/r 

^ r< l\\a \\ . i H\X7 II 

(H 2 ) ) ■ 



poo \ l/r 

y (1 + tY\\5 + d xx W(t)\\x 2 dtj < C (11^,11,1 + IIWoll^-i, 



This proves the lemma. ■ 

Lemma 5.4 For ||(0(t),W(t))||g < e, where e is sufficiently small (0 < e ^ Eq), there exists a positive 
constant C2 ^ 1 such that 

11 - •) te™) ds|lz « a '<«"-=K'))ia- (") 



Moreover, for (0 1 ,'W_ 1 ), (fl 2 , W 2 ) mi/i i/ieir norms in Z smaller than e, we have 

1 ss(t-A ( N*(gi(a),Wi(«)) -N*(0 2 ( S ),W 2 ( S )) \ 
1 NW(^ 1 ( S ),W 1 ( S ))-NW(0 2 ( S ),W 2 ( S )) 1 







(5.3) 



<c 3 (Ell&W.SM*))!!* ) IKfiiW-fiaW^W-S&W)!^ 



Proof. We start with proving the estimate (5.2). The proof is fairly straightforward. The strategy is 
to use estimates for the linear part ^M(t) in Proposition 4.2, the estimates for the nonlinear terms in 
Lemma 6.2 from the appendix, and the maximal regularity estimates in Lemma 5.1, Corollary 5.2. For 
simplicity, we denote 

N e (s) = N d (9(s), W(a)), N w (s) = N w (0(s), W(s)). 
By Lemma 6.2, we have that 

l|N*(*)||* *S —^-sWm), W(t))||| + _2_^||(e(t), W(t))|| z (l + *)||M-3£WII*», 

(1+S)2 (1 + S)4 

l|N e ( S )||, 2 *S — ^-T||(g(t), W(t))||| + _2_||(fl(t), W(t))b(l + S )||M**W(*)llx 2 , 
(1 + a)' (l + a)* 

ii^ooh* < ^-ummmiz + — ^ii(^), w(*))iu(i + «)iim«s£Wii*» 



(5.4) 



4 



|N w ( S )||x 2 *S — ^-rll(^), W(*))||| + — ^-rlKi(i), W(*))||z(l + s)\\8 + d xx W(s)\\x 2 . 

(! + «)* (1 + 8)3 
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We also exploit the linear estimates from Proposition 4.2 and obtain the following estimates. 



No 



[ Jt m {t - s)N 9 {s)ds\\ t i < C [ ||N 9 (s)||^ids, 
Jo Jo 



[ ^ i{t- s)N w {s)ds\\ i i < C [ ^^ X \ ds, 
Jo Jo (1 + t-s)* 



N 3 = (l + t)~2\\ f Jtfoo{t-s)N e {s)ds\\e°° < C(l + t)? /"* J^LWk_da, 

Jo Jo (l + t-s)2 

JV 4 = (l + t)5|| /"*^ 01 (t- S )N w (s)ds|| £ oc < C(l + t)a / l|NW(g)l|xi dg, 
Jo J 1 + t-s 

iV 5 = (l + i)f||<5 2 [' ^oo(t-s)N e (s)ds\y <C(l+i)§ f l|N (s)ll ^ 5 d5, 

Jo Jo (l + t-s)4 

iV 6 = (l + i)f ||<5 2 r^ i(t-s)N w (s)ds||^ < C(l + t)l l|NW(g)l|x2 3 dg 
Jo Jo (l + i-s)5 



i 



je w (t - s)N t '(s)ds\\ Xl < C(l + t)s /' l|N ^^. da, 
/o Jo (l + i-s)a 



JV 8 = (l + t)3|| / ^ n (t- S )N w ( s )d s || Xl <C(l + t)l / ^^feds, 
Jo J 1 + t-s 

N 9 = (l + t)i\\6+ f Jt w {t-s)n e (s)ds\\ X2 <C(l + t)f /" l|N (s)ll ^ s d g , 

Jo Jo 

iv 10 =(i+t)fii5 + r^n^-^N^dsiix^ca+t)! r ||Nw(s)||x2 3 d g , 

Jo Jo (1 + t-s)* 

= (i+*)n r^i (t- S )N e ( S )d S || Xoo <c(i+t) r&p^d Sl 

Jo Jo 1 + t - s 

iVi2 = (l + t)|| f ^ u (t- s)N™(s)ds\\ Xao ^C(l + t) T l|NW(g)l|xi T ds. 

Jo Jo (1 + t - s)(t - s)a 

At this point, we substitute (5.4) into (5.5), estimate the resulting integrals, and find 

iVj < C\\(e(t),W(t))f z , for all 1 < j < 12. 
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The calculations establishing the estimates for Ni, Nn are based on the following elementary 
integral estimates. 

1 1 C fi 1 c ri 1 

ds ^ 7- r- / ro-ds + / - — ds. 



/ (l + t- a )a(l + S )/J -(l + t )«y (l + s )/3 (1 + 7„ (1 + 5) 

For the estimate on N12, we just need to show that the following integral expression 



h(t) = (l + t) 



/ / -J- \ 

1 I.// 1 */ 1 1 X r_1 



ds + / h - rr- nr ds 



'0 (1 + t - s)(t - s)§ 1 + s Wo V( 1 + *-s)(*-s)^ + 
\ v 

has a uniform upper bound for t E (0, 00). First, for all t 6 (0, 1], there exists C > such that, 
fc(t) < 2 ^ jf (t - s)-5ds + Qf (t - s)^^ 1 ) ds^) ' J ^ C 

Second, for i G [1, 00), we have 

ft 1 1 (7 f| 1 rt 1 

(1 + i)/ TT— ds ^ r T— ds + C rd* 

Jo {l + t-s)(t-s)* 1 + 8 (l + i)l7o 1 + s 7| (l + i-a)(i-a)s 

<C|1 + / -ds ) < C. 

\ Jo (l + s)s2 J 

Similar arguments show that the second part of h(t) is also uniformly bounded on [l,oo). 

The estimates on N\, N12 bound the Z-norm of the left-hand side of (5.2), except for the maximal 
regularity component. Thus it remains to show that 

(1 + t) r \\5 + d xx W(tW X2 dt) r < C||(0(t), W(t))||!, (5.7) 



where 

W{t)= [ ^io(t-s)N e (s) + ^ n (t-s)N w (s)ds. 
Jo 

For t G [0, 1], by maximal regularity in Lemma 5.1, we have 

f (1 + ty\\5 + d xx w(t)\\ r X2 dt f 1 \\d xx w(t)\\ r X2 dt 

Jo Jo 

^^(HN^WII^ + IIN-^IUJ dt 

^C\\(9(t),mt))\\ 2 z r - 
For t E [1, 00), we split W into two parts, that is, 

-t-l r t 



(5.8) 



J + J ^ w (t-s)N e (s) + ^ 11 (t-s)lSi w ( S )ds = W 1 + W 2 . 
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By Corollary 5.2, we have 

/CO f'OQ 
(i + ty\\s + d xx w 2 (t)\\ r X2 dt <c J {i + ty\\d xx w 2 {t)\\ X2 dt 

<c jH(i + ty (\\N e (t)y + ||N w (oiu 2 ) r dt 

<c\mt) t mt))\\z- 

By similar arguments as in Lemma 5.3 and with the condition that t — s > 1, we can show that 
\\5 + d xx (uTi (< - s)N 9 (s) + ^ n (t - s)N w (s)) \\ X2 ^ - (\\N 9 (s)y + ||N w (s)||x 2 



(1 + i-s)' 

< 7 % nrll(e(<),W(*))||zf||(fl(*),W(*))||z + {1 + s)\\5 + d xx W{s)\\xX 

As a result, we obtain 

(1 + £)||5+a^i(*)IU 2 < (1 + *) J \\S+9xx {^w{t - s)N e (s) + ^ n {t - s)N w (s)) || X2 ds 
<(1 + *) f 1 " ^ nrll(0(*), W(t))|| z (||(0(t), mt))\\z + (1 + s)||M**W(*)lk 3 )d* 

./O (1 +t- + S )4 

^T-^-rlKfiW.wctjjufr, 

which immediately implies that 

(i + tyws+d^tWx^t < c|| (g(t),S£(*)) lll r - 

Together with (5.8), this establishes (5.7) and concludes the proof. 

In a completely analogous fashion, one establishes the Lipshitz estimates. ■ 
We now prove our main theorem. 

Proof of Theorem 1. The proof is a fixed-point-theorem argument. We first recall the variation of 
constant formula, 



W(t) J Wo J Jo { S) \N-(6(s),ms)), 



ds. 



Let & be the right-hand side of the formula, that is 



v w(t)y \w„J J» ' V^SW.ww), 

Assume that now the initial value is sufficiently small, that is, for some small e > 0, 

11(00; W.o)\\e 1 x{zn,T- 1 {H^)) < e. 
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If (fl(t)i W(t)) £ Z with norm smaller than e, we know that 



9_(t) 
W(t), 



e z. 



By Lemma 5.3 and 5.4, we have that 



W(t). 



W(t). 



2 



(5.9) 



Moreover, we have 



1 Iwi(t) 



2 



ll(fi,-(t),W,(t))||z HC&Ct) - fla(t),5£i(t) - Wa(t))||; 



(5.10) 

We denote -B = {V G Z | ||V||x R}, where R = mxa(2C\e,e). We now take e > small enough so 
that 2C2R < 1 and readily conclude, based on (5.9) and (5.10), that 8?{E) C B and that 8? is a strict 
contraction in B. By Banach's fixed point theorem, there is a unique fixed point of & in B, denoted as 
(9(t), W(t)). Then fflft), W(t)) is a global solution of (2.22), and if we return to the original variables, 
we obtain a global solution of (1.2) which satisfies the decay estimate in Theorem 1. This concludes 
the proof. ■ 



6 Appendix 

6.1 Estimates on nonlinear terms 

In this section, we derive the estimates on the nonlinear terms N e and N w in our normal form (2.21). 

Lemma 6.1 For || W|| y . , H^H^i < e, where e is sufficiently small(0 < e ^ 8q), there exists a nonde- 
creasing function C(e) > such that, for all 1 ^ p ^ 00, the nonlinear terms in system (2.21) have the 
following estimates. 

/ 1 \ / \ 



l N ;l E + E |0 i+fc | 3 £ KM),. 

-fc=-i Vfe=-i / Vfc=-i 

+ (E (l(*+W)i(-^)l + l(^W)i(-T)l) 

+ IIWjiiLPKtf+wjj-c-Tr)! + ywjiiip 



IN7I 



<C(e) 



E 



J+fcl 



V.Jfe=-l 



E 



+ I^II|W 



j\\LP 



+ E 



E \($+VOj + k(-«)\ + |(*+a s W) j+fc (-7r) 



0fc=-l 



+ 1 1 Wj- 1 1 ^ I (<5+ W) j- (— tt) I + ||W 2 || LP + |N e | + |N e +1 | + |Nj_!| 



(6.1) 
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Proof. We point out that throughout the proof, we repeatedly exploit the fact that the L 2 scalar 
product of an even function and an odd function are zero. We also recall that u* is even and u ad is 
odd. By equations (2. 9), (2. 10) and (2.21), we obtain 

=Ij + (II j + II I j + IV j + Vj) S*j and 

IW_( ;j 9G J \ (\rj. , T/rr. , MTTT. I TV. , 9G 3 



N ^ = \ id ~dw j ) \ VIj + VIIj + FII7j + /x > + aw**) ' where 

^ = (-1 + <W,-(x) + U j (x),n' ad (x - ; 

G j =G(%, W,-) = (W^sJ^Cs - - u ad (x))^(x - 

/,=(-^-i)(«5 + ny) j; 

//, = - (W j+1 (-7r) - Wj(— 7r), -D (u ad (vr - %) - u ad (vr))) ; 
I/Jj = (d x W j+1 (-ir) - d x Wj(-ir), Du ad (7T -9,)); 

IVj =(5 a: H i (7r) - ^H^-vr), £u ad (vr - d,)) - (H,(tt) - H^-tt), Du^tt - ^)); 

V$ =(g(%,W i + H i ),u ad (x-^)); 
VI j =A(H j -(E*6) j ); 

VIIj = - (Jllj - u[(x - 9j)0^ - (E * 9)j + (W,(x), u ad (x - 0j) - u ad (x))^(x - %)) ; 
V/JJj =(E * (<5+rw - 0)) i; 
/X, =~g{e j ,-W j + Hj) + [f'(u*0 - Oj)) - f (u*(x))] (Wj + H,); 

Xj =a(e * g)j + AWj — (e * <y + rw)j. 

We recall here that E is defined in (2.18) and point out that the term in VIIj involving Wj in fact 
cancels with a contribution from Hj. We now prove the estimate of N^. 

Estimate on I jS \Ij\ ^ C(e)\\(5 + 9)j\ + \(5-9)j\ + \\Wj\\ L p]\(8 + W)j(-iv)\. 

We first recall that Hj is defined in (2.10) and (2.13). We claim that the number Cj, appearing in the 
definition of H| as in (2.15) and (2.16), can be estimated as 



\ C j\ < C{e) 



\{8 2 9)j\ + (\(5 + 9)j\ + \(6„0)j\) J2 + \ 6 i 
^ ' k=-i 



jI||Wj||lp 



where we use notation 5 2 = 5+5-. In fact, we have 

|(0(z)(u*(z + 9j - 9 j+1 ) - u*(x + 9j - %_!)), u ad (x))| ^ C{\{5 + 9)j\ 2 + 1(5-2)/); 
|((<Kx + %) - 0(x))(u,(x + % - i+1 ) - u*(x + % - %_!)), u ad (x))| < C|%|(|(5 + ^| + \(5-9)j\); 
|(K(i + % - i+ i) + u,(x + % - %_!)), u ad (x))| ^ C|(5%|. 

We also have |Hj(x)| ^ C(e) ^\(8+6)j\ + |(5_0)j| + |%|||Wj||ipJ , from which we obtain the estimate. 

Estimate on II } : \IIj\ ^ C|6> j | 2 |(5 + W) i (-7r)|. 
This is straightforward. 

Estimate on Illji \IIIj\ < C\9j\\(5 + d x W)j(-ir)\. 
This is straightforward. 
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Estimate on IVj-. \IVj\ < C 
We first simplify IVj and obtain 

IVj =|«(tt - - uKvr - Du^tt - %)) - ^(^(tt - - u*(tt - fl^), Z?u' a d(^ - 0,)). 

Then, it is not hard to see that 



i (u*(tt - i+1 ) - u*(tt - 9j^), Du' ad (n - 9j)\ - 1 Qu^)(0f +1 - 9]^), Du'^n) 
^Cn9j +1 - 9j_ ± \ + \9j\ 2 \9j +1 - 6j-i\\, 
\ (ut&frWj-! - 9 j+1 ), -DuMej) - \ Qu*,tfe(7r)(0j +1 - 9)_ x ), Du' ad (.) 



which establishes the estimate on IVj as claimed. 

Estimate on Vji \Vj\ ^ C(e)(\(S + 9)j\ 2 + \(5-9)j\ 2 + ||W]|| LP ^. 

Noting that \Vj\ ^ C(e)||(Wj + Hj) 2 ||i P and applying the estimate of into the inequality lead to 
the above estimate. 

Estimate on S?j'. \S?j\ < C(e). 
This is straightforward. 

Combining our estimates of Ij — Vj and S^j, we obtain the first inequality in (6.1). 

Now, we have to show that the estimate of in (6.1) is true. 
Estimate on VIj\ 



\VIj\ ^ C(e) 



\(6 + 9)j\ + \(S-9)j\\ Y, \&j+k\ + \9jn\\LP 
' k=-i 



First, for / 27r-periodic and smooth, we have 

|/(z - Ox) - f(x - 9 2 ) - f'(x)(0 2 - 9x)\ < C(\9 2 - 0x\ 2 + \9 2 \\9 2 - 9x\). 
If in addition, / is odd, we have 



\m)-m)-f(Q)(e 1 -e 2 )\ ^c\9i 



11- 



The latter implies that 



\cj-\(5 2 9)j\ ^C(e)(j(«5 + 0)/ + K6-® j \ 2 + \9 j \\\W j \\ LP \ 
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Moreover, by the former inequality, we have 

\VIj\ a?(W);| 2 + 1(5-0)/ + 1^11(5+^-1 + ) + 

^C(e) 



CjAtpix - Oj) - -{6 2 0)jAiP(x) 



\(5 + 0)j\ + \(5_0)j\\ l^l + l^'lll W ill^ 
' k=-i 



Estimate on VIL: 



\VIIj\ ^ C 



E E +i%m w iii^ 



^=-i 



\k=-l 



Noting that (E * 0)j is the linear part of EL, + u*(x — Oj) — u*(x) and there is no term invovling Wj 
in VII j, we have 



\VIIj\ <C |^+i||^ +1 | + |^_i||^_i| + \0M\ + c^\x-0 3 )0j 



+ 



-(5'e)^(x)-cj^(x-d 3 ) 



where 6j = dj + (Wj(x),u ad (x - Oj) - u ad (x)). 
First, we note that 



Cj^'{x-0j)Oj ^C\0j\ 
Moreover, we claim that 



\(5 2 0)j\+{\(6 + 0)j\ + \(6„0)j\) Y, 

^ ' k=-l 



\cj\ 



\~cj - i (5%| 



\(s+o)j\ + \(s.o)j\ + ( e |0j +fc ij f E i^+fcil + i^'in w iii^ 



\k=-l 



E ' 6 'j+ fc i 

\k=~l / \k 



E i^+fei) +i^iiiw,-iu, 

:=-l / 



In fact, we have 



|<#£)(u*(* + Oj - j+1 ) - u*(x + 0, - %_!)), u ad (x))| < CI 1(5+0),- 11(5+^1 + \(6„0)j\\(6-0)j\ ) , 
(x + 0j)0j(u*{x + Oj - j+l ) - u*{x + Oj - Oj-!)), u ad (x))\ ^ C\0j\ (\{S+0)j\ + \(5.0)j 



\((0(x + Oj) - <Kx))(ii*(x + Oj ~ 6j+i) ~ + Oj - 0,-i)),u ad (x)>| < C\0j\ \{8 + 0)j\ + |(5_0),- 



\((u*(x + Oj - 6 j+1 ) + u*(x + 0j - ^_i)),u ad (x))| ^ C |(5+0),| + |(5_0) 



)3\ ' 



| <(u*(x + 0,- - J+1 ) + u*(cc + % - 0,_ 1 )),u ad (x)> + 5 2 0A ^ C 1(5+^11(5+^1 + |(5_0).,-||(5_£ 



which establishes the claim and thus the estimate on Vllj. 
Estimate on VIII,: 



\VIIIj\ ^ C[ \N°\ + |N£ +1 | + |N; 
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The calculation is straightforward using the expressions for K, and 9 a 



Estimate on IX f. 



i 

\IXj\ < C{e 



J2\(S-9) j+k \) Yl \°3+k + +1^1^-1 + |^| 2 ||W J -|| LP + |W 



|2 

J 



\k=0 / \k=-l 



The calculation is straightforward using the estimate on EL. 

dG 

Estimate on -^-Xj: 



X;| < C(e) M^l £ (|(M)i+fc| + |(5+W) J+ fc(-^)|) + |(^W j (x),u ad (x - 0j) - u ad (x)>| 



Integrating by parts, we have 

(^W j (x),U ad (x -9j)- Uad(x)) =(Wj(l),A* (u ad (x - 0y) - U ad (x))) + 

(^W i+ i(-7r) - d x Wj(-n), D(u ad (n - 9j) - u ad (vr)))- 
(W i+1 (-7r) - W^-vr), D(u' ad (7r - 6j) - u'^tt))). 

Thereofore, we have 
, 8G< 



W 3 - 



Xj\ ^ c{e)\ej\ y (l(M)j+kl + l(*+w)i+fc(-T)l) + H w ;ll^ + l(M*w)i(-70l) 



Estimate on (id — g-yy . )~ • For any £ G £°° and p G [l,oo], there exists a constant C > such that 
Combining estimates on Vlj to /-Xj, ^/ and (id — we obtain the second inequality in (6.1). 



Moreover, we have the following lemma. 

Lemma 6.2 There exist C > and rj > suc/i t/iaf, /or aZZ (0, W) G y u>i£/i iis Y -norm smaller than 
T], we have 

\\N e (s)\\ £1 < —^\\ m ,W(t))\\ 2 Y + — ^-rll(fi(t), W(t))||y(l + s)IIM**Wtollx a , 
||N*00Hp < W(t))|ft + — ^-rlKflW, W(t))|| y (l + S )||5 + ^W( S )|U 2 , 

(1 + S)2 (1 + S)2 

IIN-^IU: < T^-\\m),mt))\\ 2 Y + — ^-nrll(g(*), W(t))||y(l + S ) 11^+^ W( s ) |U 2 , 

1+S (1 + S)4 



— ^nrllCeW, W(^))lly + — ^tt II (*(*), W(*)) Ml + s)\\5 + d xx W(s)\\x 2 - 

12 



(1 + S)3 (l + *0 ! 



40 



Proof. The estimates are obtained through a direct calculation from the estimates in Lemma 6.1. We 
sketch the computation for ||N e (s)||^i, and the others follow similarly. 

First, for terms only involving 9, we notice that 

£l(M)jl 2 = -E^% ^ \\d\\ i2 \\5 2 ey < —J-^i^wct))^, 

96Z igZ (1 + Sj 2 



E>il 3 l(M)il < l|0||2»||£IHIMI* < ||0||^||0||f||<5 2 0||| < ?T ^||(0(t),W(*))llr- 
iez 11 + sj 

Second, for terms involving ^W, we observe that 



E |eil|(M«w)i(-T)l < E ( T (^ w i+i(^) - ^(i)) dx) 



2 



< V27r\\0\\ e 2\\5 + d xx W\\x 2 



< -^-s\\m), W(t))||y(l + S )||M-W(*)lk- 

(1 + S )4 

Similarly, for X)jez |^j||(^-i-W)j(— 7r)|, we have 

Sl^ll(*+W)i(-7r)l < ^\\0y\\6+d x W\\x*- 

jez 

Using the "homogeneous matching boundary conditions " (2.11), we have 



\\5 + d x W\\x 2 = ( - E r ( S +^0j ( x ) ( J +^W) . (x)dx 



i 



< II^+wii^jm^wii^ 

< \\S+W\\x 2 + \\5 + d xx W\\x 2 . 
We plug the latter estimate into the former one and obtain that 

ENI(<*+wm-7t)| < V2^\\ey ( \\6+w\\x a + \\s+d xx w\\x 2 ) 



< 7^%-||(*(*), W(*))IIy + T^rlKgW, W(t))||y(l + s )||« + flUW( s )||* a . 

(1 + S)5 (1 + S)4 

For Sjez ||Lp|(5+W)j(— 7r)|, we take p = 2 and follow steps as above, obtaining the following 
estimate. 



E llw.-lbKMVM-Tr)! < -^-ik^), w(t))|| y +-^U|(0(i), S£(t))||y(i+*)||« + Q«5£(*)l|j6. 

jez + + 

For Sjez we take p = 1 and obtain that 

jez (1 + s) 2 

Combining the above estimate, we establish the first inequality in the lemma. ■ 
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6.2 Bloch wave decomposition 

In this section, we present the Bloch wave decomposition of the linear operator A. We first recall that 
A, as in (3.1), is defined as 

A: (H 2 (R)) n — ► (L 2 (M)) n 

v i — > Dd xx v - f'(u*)v. 

We introduce the direct integral [17, XIII. 16.] 

39: L 2 (Ti, (L 2 (T 27r )) n ) — ► (tf(R))« 

U(a,x) J CT6Ti e--U(a,-)da- 1 • ^ 

The direct interal is an isometric isomorphism with inverse 

m~ x : (L 2 (R)) n — ► L 2 (T 1; (L 2 (T 27r ))") 
u(z) ■— > ^£ teZm e i£ - x u(cr + f). 

The following result from [20, 14] characterizes the Bloch wave decomposition of A. 

Theorem 2 (Bloch wave decomposition) The linear operator A is diagonal in Bloch wave space. 
To be precise, 



X ASS = A= / B(a)da, (6.3) 



where by A = J 2 1 B(o~)da, we mean that, given any u G L 2 (Ti, (L 2 (T2,r)) n ), 

— 2 

(Au)(a) = B{a)u(a), a.e. a G [--, -]. 
Moreover, we have the following spectral mapping property. 

spec(A) = spec(A) = |J spec(B(a)). (6.4) 

6.3 Spectral properties of {A ch (a)} cre ^_i ij 

We recall that ^4 c h(u) is defined in (3.9) as A c h(a) = ^ n B(a)^'~ 1 and Y q in (3.15) for 1 ^ q ^ oo. 
We are concerned with their spectral properties as unbounded operators in Y q , which is useful for the 
derivation of the estimates for M(t,o~) as defined in (3.23). 

We first show the well-definedness of A c h(a) in Y q in the following lemma. 

Lemma 6.3 For any given a G [— I, \], A^a) is an unbounded closed operator in Y 2 , that is, 

A ch (a): 3> 2 (A ch (o-))cY 2 — > Y 2 

w i — > {-(a + l) 2 Dwi + Y,kez h £-k w k}eez, 

where & 2 (A ch (a)) = {w G Y 2 | {(1 + m 2 )w m } meZ G Y 2 } and h = {h £ } te z = ± JZj' (u*(x))e- ikx dx. 
Moreover, A^a) can naturally be considered as an unbounded closed operator inY q , with ^(A^a)) = 
{w G Y q | {(1 + m 2 )w m } meZ G Y q }, for all 1 ^ q < oo. 
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Proof. The expression for A c ^(a) in Y<± follows from a direct calculation. The extension to Y„ fol- 
lows from the fact that the set {w £ | w has finitely many nonzero entries} is dense in Y q and 
%{A ch {a)), for all q £ [1, oo]. ■ 

We then have the following proposition. 

Proposition 6.4 For any fixed a £ [— i, i] and j> £ [l,oo], A c h(cr) defined in Y q is sectorial and has 
compact resolvent. In fact, there exist C > 0, ui £ (tt/2,tt) and Ao £ R, independent of a and q, such 
that the sector S(\q,u)) = {A £ C | ^ | arg(A — Ao) | ^w,A/ Ao} Q p{A c \ l {a)) and 

|||(l ch (a) - A) -1 |||y 9 < C|A - A r\ for all A £ S(X ,uj),a £ [-1 1] and q £ [l,oo]. (6.6) 

Moreover, for any fixed a £ [— o> o]> ^ e s P ec ~k rum °f A c h(c) independent of the choice of its under- 
lying space Y q and thus denoted as spec(A c h(cr)), for any q £ [l,oo], with spec(A c h(cx)) = spec(B(cr)) 
consisting only of isolated eigenvalues with finite multiplicity. 

Proof. We view A c h(<r) as a perturbation of the Laplacian in the discrete Fourier space, that is, 

A ch (a) =L(a) + H, 

where L(a)w_ = {—(a + Q 2 Dwf}fa and Hw = {X^/tgz h^-fc w fc}tez- It is straightforward to verify 
that the proposition holds for the Laplacian L{a). We only have to show that the perturbation H is 
good enough to preserve these properties. Noting that H £ J£ ((£ p ) n ) for any p £ [1, oo] with its norm 
uniformly bounded, we have, for A £ p(L(a)), |A| sufficiently large, 

(A ch (a) - A)" 1 = (L(a) + H- A)" 1 = (L(a) - A)' 1 (id + H(L(a) - A) -1 ) -1 . (6.7) 

All assertions in the proposition easily follows from this expression (6.7), except for the fact that the 
spectrum of A c ^(a) is independent of q. 

To prove this property, we denote the spectrum of A c h(cr) defined on Y q as spec(A c h(<r), q), which 
consists of eigenvalues with finite multiplicity, accumulating at infinity, only. Given any eigenfunction 
v = {vj}j 6 z, v belongs to Dgerioo]^? since v are smooth, that is, Vj decays algebraically with any 
rate. This establishes spec(A(<r), g) = spec(A(a) , p) , for any p, q £ [l,oo]. ■ 

6.4 Perturbation results 

We apply perturbation theory to the Bloch wave operator B(o~) for a near and obtain more detailed 
spectral information, including the Taylor expansion of d in Hypotheses 1.2. 

To this end, we define 

F : [-i,|]xCxiJ| — > I? 

(a,X,w) — > (S( fT )-A)(w + u'J, 

where Hj_ = {w £ (H 2 (T2 7T )) n | (w,u^) = 0}. A standard implicit-function-theorem argument shows 
that there are a small neighborhood of a at the origin and a smooth function (A(cr), w(<r)) with 
(A(cx), w(<r)) = on this neighborhood such that F(a, X(a),w(a)) = 0. We denote e(cr) = u^. + w(cr). 
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Similarly, replacing B{a) with its adjoint B*(a), we obtain a smooth continuation of u a( j, denoted 
as e*(c). Without loss of generality, we can assume that (e(a"), e* (cr)) = 1. Moreover, we have the 
following proposition. 

Proposition 6.5 There exist positive numbers 70 and 71 such that for any \cr\ ^ 70 in R, B(o) has 
only one simple eigenvalue within the strip | Re A| ^ 71 in C, which is exactly the continuation A(cr) of 
the eigenvalue A(0) = 0. Moreover, A(cr) has the Taylor expansion, 

A(o-) = -da 2 + 0(H 3 ), 

where — 7i/4 ^ —2da 2 < ReX(a) < — \<J 2 , for all a G [— 70,70] «^ 

, <9 2 e(0, x) , . . „ ... 
d = -<2i - <(x), Z?u ad (x)>. 

Proof. We first derive the explicit expression of d. To do that, taking first and second derivative with 
respect to a of F(a, A(<r), w(<r)) = 0, taking the inner product of the derivatives with u a( j and letting 
a = 0, we have 

A'(0) = (B(0)d a e(0,x) +2iD<(aO,u ad (aO>, 

A"(0) = (5(0)^e(0,x) + (4i^^ - 2A'(0))5 CT e(0, x) - 2Z><(x),u ad (x)). 

Noting that span{u a( j} _L Rg(B(0)) and the inner product of an even function and an odd function is 
always 0, we have 

A'(0) = 0, A"(0) = 2(2i ^ ( ^ x) - u;(x),£> Uad (x)). 

It remains to prove the uniqueness of the eigenvalue of B{a) in a vertical strip centered at the origin 
for sufficiently small a. First, there is no eigenvalue within the strip far away from the origin due to 
the fact that, by Proposition 6.4, spec(B(a)) is in the same sector for every a G [— |,o]- Secondly, 
the uniqueness within a small neighborhood of the origin follows from the above perturbation results. 
For the region inbetween, compactness and the local robustness of resolvent guarantee the absence of 
eigenvalues within this area. ■ 

Remark 6.6 (i) We stress that we may choose 70 as small as desired. 

(ii) The uniqueness implies that, for \a\ sufficiently small, A(cr) is a real number since its complex 
conjugate is also an eigenvalue. 



6.5 Properties of analytic semigroups ^e Ach ^ t } a 
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In this section, we will derive various estimates on e jich ^ a ' t . We first note that by [5, 1.4] the interpolation 
space S) q {A c \ l {a) a ) is independent of a, 



%(A ch (a) a ) = {w G Y q I {(1 + m 2 ) a w m } mez G Y q } =: Y q a , ||w||y- = ||{(1 + m 2 ) a w m } m& \\ Yq . 

We then recall the definitions of Y q>c (a), Y q>B (a), A c (a) and A s (a) from (3.18). We now have the 
following proposition. 
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Proposition 6.7 For every q G [1, +00] and a > 0, there exist positive constants e 6 (0,1), 72, C(q), 
C(a) and C(a,q) such that 

e-% aH , for all \a\ <7o,i ^ 0, 
C{a)e~^ a2t , for all \a\ ^ j ,t ^ 0, 
C(q)e~~^\ for all \a\ < 70, t ^ 0, 
C(a,g)t~ a e~ 7lt/2 , /or oM |ct| < 7o ,t > 0, 
C{q)Q- eda2 \ for all \a\ < 70, i > 0, 
C(9)e" 72 *, /or a// 7o < |ct| < -,t ^ 0, 
C(a,p)t- Q e- 72 *, for all 7o < |a| < -,t > 0. 

Proof. We first derive estimates for the case \a\ ^ 70- For A c (a), we have e Ac{ - a)t = e A(<j)t . The first 
two inequalities follow directly from the fact that ReA(a) < — |<7 2 and e(cr) is smooth, by Proposition 
6.5, for \o~\ ^ 7o- 

For A s (a), by Proposition 6.4 and 6.5, for any a G (—70,70) an d q G [l,oo], 

^ 7l 7T 

spec(A s (cr),g) C C\5( — — ,uj), where u) G ( — , 7r). 
Moreover, for every q G [1, +00], there exists a positive constant C(q) such that 

\UMo-) - A) -1 III Vg, s (or) < C(q)\X + ^r 1 , for all \a\ ^ 70 and A G S(- 

Thus, by [5, Thin. 1.3.4, 1.4.3], we immediately obtain the two inequalities for A s (o~). The first inequality 
on A c h(a) follows directly by combining the first inequality for A c (a) and the first inequality for A s (a). 

We now derive the estimates for the case 70 < |o~| ^ \. By a similar analysis as in Proposition 6.5, 
there exists a positive constant 72 such that 

Re(specA c h(o")) < —272, for all 70 < \a\ < -. 
It is then not hard to conclude that 

^ 7T 

spec(A ch (o-)) C Cy> (-72,0)1), where cji G ( — , 7r). 

Moreover, for every q G [1, +00], there exists a positive constant C(q) such that 

|||(l ch (o-) - A)" 1 !!!^ < C(g)|A + 7 2 r 1 , for all 7o < |cr| < - and A G £(-72,0)1). 

Therefore, again by [5, Thm.1.3.4, 1.4.3], we immediately obtain the last two inequalities for A,±(a), 
which concludes the proof. ■ 



|pA:(o-)t||| / < 



9,c(ct) 



le As ( ff )*lll , ^ < 



3 ^chW*||| < 



^A ch (a)t\ 



\Y^Y« 



45 



References 

[1] J. Bricmont and A. Kupiainen, Renormalization group and the Ginzburg-Landau equation, 
Comm. Math. Phys., 150 (1992), pp. 193-208. 

[2] , Stability of moving fronts in the Ginzburg-Landau equation, Comm. Math. Phys., 159 (1994), 

pp. 287-318. 

[3] K. Deng and H. A. Levine, The role of critical exponents in blow-up theorems: the sequel, J. 
Math. Anal. Appl., 243 (2000), pp. 85-126. 

[4] T. Gallay and A. Scheel, Diffusive stability of oscillations in reaction- diffusion systems, Trans. 
Amer. Math. Soc, 363 (2011), pp. 2571-2598. 

[5] D. Henry, Geometric theory of semilinear parabolic equations, vol. 840 of Lecture Notes in Math- 
ematics, Springer-Verlag, Berlin, 1981. 

[6] M. A. Herrero and J. J. L. Velazquez, Some results on blow up for semilinear parabolic 
problems, in Degenerate diffusions (Minneapolis, MN, 1991), vol. 47 of IMA Vol. Math. Appl., 
Springer, New York, 1993, pp. 105-125. 

[7] M. A. Johnson, Nonlinear stability of periodic traveling wave solutions of the generalized 
Korteweg-de Vries equation, SIAM J. Math. Anal, 41 (2009), pp. 1921-1947. 

[8] M. A. Johnson and K. Zumbrun, Nonlinear stability of periodic traveling wave solutions of 
systems of viscous conservation laws in the generic case, J. Differential Equations, 249 (2010), 
pp. 1213-1240. 

[9] , Nonlinear stability of periodic traveling-wave solutions of viscous conservation laws in di- 
mensions one and two, SIAM J. Appl. Dyn. Syst., 10 (2011), pp. 189-211. 

[10] , Nonlinear stability of spatially-periodic traveling-wave solutions of systems of reaction- 
diffusion equations, Ann. Inst. H. Poincare Anal. Non Lineaire, 28 (2011), pp. 471-483. 

[11] M. A. Johnson, K. Zumbrun, and P. Noble, Nonlinear stability of viscous roll waves, SIAM 
J. Math. Anal., 43 (2011), pp. 577-611. 

[12] A. Lunardi, Analytic semigroups and optimal regularity in parabolic problems, Progress in Non- 
linear Differential Equations and their Applications, 16, Birkhauser Verlag, Basel, 1995. 

[13] A. MlELKE, Uber maximale L p -Regularitat fur Differentialgleichungen in Banach- und Hilbert- 
Raumen, Math. Ann., 277 (1987), pp. 121-133. 

[14] , Instability and stability of rolls in the Swift- Hohenberg equation, Comm. Math. Phys., 189 

(1997), pp. 829-853. 

[15] J. D. Murray, Mathematical biology. I, vol. 17 of Interdisciplinary Applied Mathematics, 
Springer-Verlag, New York, third ed., 2002. An introduction. 



46 



[16] , Mathematical biology. II, vol. 18 of Interdisciplinary Applied Mathematics, Springer- Verlag, 

New York, third ed., 2003. Spatial models and biomedical applications. 

[17] M. Reed and B. Simon, Methods of modern mathematical physics. IV. Analysis of operators, 
Academic Press [Harcourt Brace Jovanovich Publishers], New York, 1978. 

[18] B. Sandstede and A. Scheel, On the stability of periodic travelling waves with large spatial 
period, J. Differential Equations, 172 (2001), pp. 134-188. 

[19] B. Sandstede, A. Scheel, G. Schneider, and H. Uecker, Diffusive mixing of periodic wave 
trains in reaction- diffusion systems, J. Differential Equations, 252 (2012), pp. 3541-3574. 

[20] B. SCARPELLINI, P -perturbations of periodic equilibria of reaction diffusion systems, NoDEA Non- 
linear Differential Equations Appl., 1 (1994), pp. 281-311. 

[21] G. Schneider, Diffusive stability of spatial periodic solutions of the Swift- Hohenberg equation, 
Comm. Math. Phys., 178 (1996), pp. 679-702. 

[22] A. Turing, The chemical basis of morphogenesis, Philosophical Transactions of the Royal Society 
of London. Series B, Biological Sciences, 237 (1952), pp. 37-72. 

[23] H. Uecker, Diffusive stability of rolls in the two-dimensional real and complex Swift- Hohenberg 
equation, Comm. Partial Differential Equations, 24 (1999), pp. 2109-2146. 

[24] S. Zelik and A. Mielke, Multi-pulse evolution and space-time chaos in dissipative systems, 
Mem. Amer. Math. Soc, 198 (2009), pp. vi+97. 



47 



